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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement 
be  converted  to  metric  (SI)  units  as 


Multiply 

degrees 

feet 

inches 


0.01745329 

0.3048 

25.4 


used  in  this  report  can 
follows : 

To  Obtain 
radians 
metres 
millimetres 
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RESONANT  COLUMN  TEST 


1.  INTRODUCTION 

T he st  methods  cover  the  determination  of  the  shear  modulus, 
shear  damping,  rod  modulus  (commonly  referred  to  as  Koung's 
modulus),  and  rod  damping  of  cylindrical  specimens  of  soil  in  the 
undisturbed  and  remolded  conditions  by  vibration  using  the 
re  sonant  column.  The  vibration  of  the  specimen  may  be  superposed 
on  a controlled  ambient  state  of  stress  in  tne  specimen.  The 
vibration  apparatus  and  specimen  may  be  enclosed  in  a trinxial 
chamber  and  subjected  to  an  all  around  pressure  and  axial  load. 
In  addition,  the  specimen  may  be  subjected  to  other  controlled 
environment  conditions  (for  example  temperature).  These  methods 
for  modulus  and  damping  determination  are  considered 
nondestructive  when  the  strain  amplitudes  of  vibration  are  less 
tnan  10f-4  in. /in.  and  many  measurements  may  be  made  on  the  same 
specimen  at  a given  ambient  stress  and  at  various  other  states  of 
ambient  stress.  These  methods  cover  the  determination  of  the 
modulus  and  damping  , the  necessary  vibration,  and  specimen 
preparation  procedures  related  to  the  vibration,  etc.,  and  the 
application,  measurement  and  control  of  the  ambient  stress.  The 
application,  measurement,  and  control  of  ambient  stress  is  quite 
similar  to  that  given  in  Appendix  X:  Triaxial  Compression  Tests? 
The  description  of  apparatus,  apparatus  calibration,  and 
calculations  (including  the  enclosed  computer  program)  have  been 

* Engineering  Manual  1110-2-1906  , 'Laboratory  Soils  Testing," 
November  1970. 
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adapted  from  Drnevich  et  al.* 


2.  SIGNIFICANCE 

The  modulus  ana  damping  of  a given  soil,  as  measured  by  the 
resonant  column,  depend  upon  the  strain  amplitude  ot  vibration, 
the  ambient  state  of  effective  stress,  and  *he  void  ratio  of  the 
soil  as  well  as  other,  less  significant  factors.  The 
applicability  of  the  results  to  the  rield  conditions  will  depend 
on  the  degree  to  which  the  application  and  control  of  the  ambient 
stresses  and  the  void  ratio,  as  well  as  other  parameters, 
duplicate  field  conditions. 

3.  DEFINITIONS 

a.  Resonant  Column.  - A cylindrical  specimen  or  column  of 
soil  that  has  platens  attached  to  each  end  as  shown  in  Figure  1, 
A sinusoidal  vibration  excitation  device  is  attached  to  the 
active  end  platen.  Tne  other  end  is  the  passive  end  platen.  It 
may  be  rigidly  fixed  or  its  mass  and  rotational  inertia  must  be 
known.  The  vibration  excitation  device  may  incorporate  springs 
and  dashpots  of  known  characteristics  connected  to  the  active  end 
platen.  Vibration  excitation  may  be  longitudinal  or  torsional. 
A given  apparatus  may  have  the  capability  of  one  or  the  other,  or 
both.  The  mass  and  rotational  inertia  of  the  active  end  platen 
and  portions  of  the  vibration  excitation  device  moving  with  it 


* Drnevich,  V.F.,  Hardin,  B.O.,  and  Shippy,  D. J. , "Modulus  and 
Damping  of  Soils  by  the  Resonant  Column  Method,"  ASTM  Symposium 
on  Dynamic  Soil  and  Rock  Testing  in  the  Field  and  Laboratory  for 
Seismic  Studies,  Denver,  Colorado,  June,  1H77. 
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(2)  For  the  case  where  platen  mass  (or  passive  end  platen 
rotational  inertia)  is  greater  than  100  times  the  corresponding 
value  of  the  specimen  and  is  not  rigidly  fixed,  resonance  is  the 
lowest  frequency  for  which  the  sinusoidal  excitation  force  (or 
moment)  is  100  degrees  out  of  phase  with  the  velocity  of  the 
active  end  platen. 
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b.  System  Kesonant  Frequency.  - The  definition  of 
resonance  depends  on  both  apparatus  and  specimen  character 
(1)  Where  the  passive  end  platen  is  fixed,  motion 
active  end  is  used  to  establish  resonance  which  is  defined 
lowest  frequency  tor  which  the  sinusoidal  excitation  fo 
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c.  Ambient  Stress.  - Static  confining  stresses  (whether 
isotropic  or  anisotropic'  applied  to  the  specimen  during  the  test 
are  termed  ambient  stresses.  The  ambient  stress  at  the  time  of 
measurement  of  the  system  resonant  frequency  and  system  damping 
shall  be  measured  and  recorded  in  accordance  with  Section  7. 

d.  Moduli  and  Damping  Capacities.  - (1)  Young's  modulus 
(herein  called  rod  modulus),  E,  is  determined  from  longitudinal 
vibration  ana  the  shear  modulus,  G,  is  determined  from  torsional 
vibration.  The  rod  and  shear  moduli  shall  be  defined  as  the 
elastic  moduli  oi  a unirorm,  linearly  viscoelastic  specimen  of 
the  same  mass  density  and  dimensions  as  the  soil  specimen, 
necessary  to  produce  a resonant  column  having  the  measured  system 
resonant  frequency.  The  stress-strain  relation  for  a steady- 
state  vibration  in  the  resonant  column  is  a hysteresis  loop. 
These  moduli  will  correspond  to  the  slope  of  a line  through  the 
end  points  of  the  hysteresis  loop.  Section  6 provides  for 
computation  of  shear  and  rod  moduli  from  the  measured  system 
longitudinal  and  torsional  resonant  frequencies,  respectively. 

(2)  The  energy  dissipated  by  the  system  is  a measure  of  the 
damping  capacity  of  the  soil.  Damping  will  be  defined  by  the  rod 
damping  ratio,  DL,  ana  the  shear  damping  ratio,  D T,  which  are 
analogous  to  the  critical  viscous  damping  ratio,  c/ccr,  for  a 
single  degree  of  freedom  system.  The  damping  ratios  for  rod 
compression  shall  be  defined  by 

Dl  = 0 . 5 ( i|uj/E) 
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where:  i|  = viscous  coefficient  for  rod  compression 

uu  = circular  resonant  frequency 

and  by 

D = 0.5(nuu/G) 

where:  = viscous  coefficient  for  shear. 

Values  for  damping  determined  in  this  way  will  correspond  to  the 
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it  is  common 

to  use  complex 

moduli 

to 

ex  press 

both 

modulus 

and  damping. 

The  complex  rod 

modu lus 

is 

given  by 

E*  = E ( 1 t 2iD  ) 

where:  i = square  root  of  -1, 

and  the  complex  shear  modulus  is  given  by 

G*  = G ( 1 ♦ 2iD  ) 

e.  Cyclic  Strain.  - (1)  The  sinusoidal  vibration  excitation 
causes  sinusoidal  vibratory  motion  having  equal  positive  and 
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negative  magnitudes. 

(2)  For  longitudinal  excitation,  the  strain,  e,  for  each 
half  cycle  (positive  half  or  negative  half)  is  the  average  axial 
strain  in  the  entire  specimen. 

(3)  For  torsional  excitation,  the  strain,  Y,  for  each  half 
cycle  is  the  average  shear  strain  in  the  specimen.  In  the  case 
of  torsion,  shear  strain  in  each  cross  section  varies  from  zero 
along  the  axis  of  rotation  to  a maximum  at  the  perimeter  of  the 
specimen  and  the  average  shear  strain  for  each  cross  section 
occurs  at  a radius  egual  to  2/3  the  radius  of  the  specimen. 

f.  Apparatus  Model  and  Constants.  - (1)  The  rigidity  and  mass 
distribution  of  the  resonant  column  shall  be  as  required  in 
Section  5 in  order  for  the  resonant  column  system  to  be 
accurately -represented  by  the  model  shown  in  Figure  1. 

(2)  The  apparatus  constants  are:  the  mass  of  the  passive 
end  platen,  M ^ including  the  mass  of  all  attachments  rigidly 
connected  to  it;  the  rotational  inertia  of  the  passive  end 
platen,  Jp,  including  the  rotational  inertia  of  all  attachments 


rigidly 

connected 

to 

it 

; sim 

ilar 

mass,  N . , and  rotational 
A 

inertia. 

JA,  for 

the 

ac 

ti  ve 

end 

platen  and  all  attachments 

rigidly 

connected 

to 

it 

such 

as 

portions  of  the  vibration 

excitation  device;  the  spring  and  damping  constants  for  both  the 
longitudinal  and  torsional  springs  and  dashpots  (K ^ , K ^ , ADC^, 
ADC  ) ; the  apparatus  resonant  frequencies  for  longitudinal 
vibration,  f0L * and  torsional  vioration,  f0T  5 the  force/current 
constant,  FCF,  relating  applied  vibratory  force  to  the  current 
applied  to  the  longitudinal  excitation  device;  the  torque/current 
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constant,  TCF,  relating  vibratory  torque  to  the  current  applied 
to  the  torsional  excitation  device;  and  factors  (LCF^ , RCF^  , LC^. 
, RCFp  ) relating  the  transducer  outputs  to  active  and  passive 
end,  longitudinal  and  rotational  motion. 

4.  APPARATUS 

a.  General.  - The  complete  test  apparatus  includes:  the 
platens  for  holding  the  specimen  in  the  pressure  cell,  the 
vibration  excitation  device,  transducers  for  measuring  the 
response,  the  control  and  readout  instrumentation,  and  auxiliary 
equipment  for  specimen  preparation. 

b.  Specimen  Platens.  - (1)  Both  the  active  end  and  passive 
end  platens  shall  be  constructed  of  noncorrosive  material  having 
a modulus  at  least  ten  times  the  modulus  of  the  material  to  be 
tested.  Each  platen  shall  have  a circular  cross  section  and  a 
plane  surface  of  contact  with  the  specimen,  except  that  the  plane 
surface  of  contact  may  be  roughened  to  provide  for  more  efficient 
coupling  with  the  ends  of  the  specimen.  Porous  stones  or  porous 
metallic  discs  of  non-cor rosive  materials  may  be  used  for  this 
purpose. 

(2)  Pore  water  drainage  lines  may  be  connected  to  one  or 
both  platens;  however,  it  is  important  that  the  lines,  when 
pressurized,  do  not  affect  apparatus  resonant  frequency  and 
damping  or  that  their  effect  be  accounted  for  in  calibration. 

(3)  Pore  pressure  transducers  may  be  mounted  directly  in 
the  platens  or  may  be  mounted  in  the  pore  water  line(s).  When 
drained  tests  are  performed,  the  transducer  (s)  measure  the 
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ambient  back  pressures  if  back  pressures  are  used.  fn  undrained 


tests  where  pore  water  drainage  valves  are  closed,  the 
transducer  (s)  measure  the  ambient  pore  pressure (s)  at  the 
specimen  boundaries.  The  transducer  (s)  cannot  accurately  measure 
dynamic  pore  pressures  unless  special  precautions  are  taken  in 
the  design  of  the  pore  pressure  measuring  system  and  the  pore 
pressure  equalization  time  for  the  specimen  is  less  than  0.01 
times  the  period  of  the  applied  vibration  (inverse  of  the  system 
resonant  frequency) . For  details  of  accurate  pore  pressure 
measurement  see  Bishop  and  Henkel.* 

(4)  The  diameter  of  the  platens  shall  be  equal  to  or 
greater  than  the  diameter  of  the  specimen.  The  active  end  platen 
may  have  a portion  of  the  excitation  device,  springs,  and 
dashpots  connected  to  it.  These  must  be  rigidly  connected  and 
are  to  be  considered  part  of  the  platen.  The  passive  end  platen 
may  have  a mass  rigidly  attached  to  it  or  it  may  be  rigidly 
fixed.  The  construction  of  the  platens  shall  be  such  that  they 
are  essentially  rigid  with  respect  to  the  specimen. 

c.  Vibration  Excitation  Device.  - (1)  This  shall  be  an 
electromagnetic  device  capable  of  applying  a sinusoidal 
longitudinal  vibration  and/or  torsional  vibration  to  the  active 
end  platen  to  which  it  is  rigidly  coupled. 

(2)  The  frequency  of  excitation  shall  be  adjustable  and 
controlled  to  within  0.  1 percent.  The  rigidity  and  mass 
distribution  of  the  vibration  excitation  device  when  fastened  to 


* Bishop,  A.H.,  and  Henkel,  D.J.,"The  Triaxial  Test,"  Edward 
Arnold  Publishers,  London,  1962,  228  pages. 
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the  active  end  platen  shall  be  such  that  it  can  be  accurately 

represented  as  a rigid  mass  that  may  be  attached  to  weightless 

springs  and  dashpots  as  shown  in  Figure  1.  If  weightless  springs 

are  used,  the  excitation  device  and  active  end  platen  (without 

the  specimen  in  place)  form  a two  deqree-of-f reedom  system  (one- 

degree-of-f reedom  system  for  devices  designed  for  only 

longitudinal  or  only  torsional  motion) , having  undamped  natural 

frequencies  for  longitudinal  motion,  foL  , and  torsional  motion, 

f . The  device  shall  be  constructed  such  that  these  modes  of 
oT 

vibration  are  uncoupled. 

(3)  The  excitation  device  shall  have  a means  of  measuring 
the  applied  current  that  has  at  least  a 5 percent  accuracy.  The 
voltage  drop  across  a fixed,  temperature  and  frequency  stable, 
power  resistor  in  series  with  the  excitation  device  may  be  used 
for  this  purpose.  The  force/current  and  torque/current  factors 
for  the  vibration  excitation  devices  must  be  linear  within  5 
percent  for  the  entire  range  of  operating  frequencies  anticipated 
when  testing  soils. 

d.  Sine  Wa ve  Generator.  - (1)  An  electric  instrument  capable 

of  producing  a sinusoidal  current  with  a means  of  adjusting  the 
frequency  over  the  entire  range  of  operating  frequencies 
anticipated.  This  instrument  shall  provide  sufficient  power  to 
produce  the  required  vibration  amplitude,  or  its  output  may  be 
electronically  amplified  to  provide  sufficient  power.  The  total 
distortion  of  the  signal  applied  to  the  excitation  device  shall 
be  less  thar  3 percent.  The  sensitivity  of  the  trejuency 
adjustment  shall  be  less  than  0.1  percent. 
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e.  Vibration  Measuring  Devices  and  Readout  Instruments.  - (1) 
The  vibration  measuring  devices  shall  be  acceleration,  velocity, 
or  displacement  transducers  that  can  accurately  measure  the 
motion  of  the  active  and  passive  end  platens.  A transducer  is 
not  needed  on  the  passive  end  platen  if  it  rigidly  fixed. 

(2)  On  each  platen,  one  transducer  shall  produce  a 
calibrated  electrical  output  that  is  proportional  to  the 
longitudinal  acceleration,  velocity,  or  displacement  of  that 
platen  (not  reguired  for  torsion  only  apparatus) . The  other 
transducer  shall  produce  a calibrated  electrical  output  that  is 
proportional  to  the  rotational  acceleration,  velocity  or 
displacement  (not  reguired  for  longitudinal  only  apparatus).  The 
readout  instrument  and  transducers  shall  have  a sensitivity  such 
that  a displacement  of  10f-5  in.  (2.b  x 10|-7  n)  and  a rotation 
of  10f-5  radian  can  be  measured  with  10  percent  accuracy  for  the 
entire  range  of  frequency  anticipated. 

(3)  It  is  also  necessary  to  have  an  x-y-time  oscilloscope 
available  for  observing  signal  waveforms  and  for  establishing  the 
system  resonant  frequency.  This  oscilloscope  must  have  at  least 
one  amplifier  (vertical  or  horizontal)  with  sufficient  gain  to 
observe  the  motion  transducer  output  over  the  entire  range  of 
output  voltages  and  frequencies  anticipated. 

(4)  For  measurement  of  damping  by  the  tree  vibration 
method,  and  for  calibration  of  the  apparatus  damping,  the  readout 
instrument  shall  be  capable  of  recording  the  decay  of  free 
vibration.  Either  a strip  chart  recorder  with  appropriate 
response  time  and  chart  speed  or  an  oscilloscope  and  camera  may 
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be  used  for  this  purpose. 

f.  Support  for  Vibration  Excitation  Device.  - For  the  special 
case  where  the  vibration  excitation  device  and  active  end  platen 
are  placed  on  top  of  the  specimen,  it  may  be  necessary  to  support 
all  or  a portion  of  the  weight  of  the  platen  and  excitation 
device  to  prevent  excessive  axial  stress  or  compressive  failure 
of  the  specimen.  This  support  may  be  provided  by  a spring, 
counter  balance  weignts  or  pneumatic  cylinder  as  long  as  the 
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described  in  Appendices  X or  XI  may  be  used  for  specimen 
preparation  and  application  of  ambient  stresses-,  Additional 
miscellaneous  apparatus  and  tools  may  be  used  for  these  purposes 
as  required. 

5.  APPARATUS  CALIBRATION 

a.  Notion  Transducers.  - (1)  Notion  transducers  shall  be 
calibrated  with  each  other  and  with  an  independent  method  to 
insure  calibration  accuracy  within  5 percent. 

(2)  Linear  motion  transducers  whose  axes  are  located  fixed 
distances  from  the  axis  of  rotation  may  be  used  to  measure 
rotational  motion  if  the  cross-axis  sensitivities  of  the 
transducers  are  less  than  5 percent.  For  this  case  the  distance 
between  the  axis  of  rotation  and  the  transducer  axes  shall  he 
known  to  within  5 percent. 

(3)  The  calibration  factors  for  longitudinal  motion  shall 
be  expressed  in  terms  of  peak  meters/peak  volt.  The  calibration 
factors  for  rotational  motion  shall  be  expressed  in  terms  of  peak 
radians/peak  volt.  This  means  that  for  velocity  and  acceleration 
transducers,  the  vibration  frequency  shall  be  included  as  a term 
in  *he  calibration  factor.  For  velocity  transducers,  the 
calibration  factors  are  given  by 

Displ.  Cal.  Fact.  = Vel.  Cal.  Fact./(2wf) 

and  for  acceleration  transducers,  the  calibration  factors  are 
given  by 
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Displ.  Cal.  Fact 


Accel.  Cal.  Fact ./ (2nt ) | 2 


Thus,  for  velocity  and  acceleration  transducers,  the  calibration 
factors  will  not  be  constants  but  will  vary  with  measured 
frequency,  f. 

(4)  Calibration  factors  tor  longitudinal  motion  are  given 
by  the  symbol  LCF  with  a subscript  A or  P denoting  whether  the 
transducer  is  associated  with  the  active  end  platen  or  passive 
end  platen.  Likewise,  the  calibration  factors  for  rotational 
motion  will  be  given  by  the  symbol  HCF  and  will  have  subscripts  A 
or  P depending  on  their  location. 

b.  Passive  F.nd  Platen  (lass  and  Rotational  Inertia.  - (1)  The 
mass  and  rotational  inertia  of  the  passive  end  platen  shall  be 
determined  with  all  transducers  and  other  rigid  attachments 
securely  in  place.  The  mass,  Sp,  is  calculated  from 


where:  Wp  = the  weight  of  the  passive  end  platen  and  its 
attachments 

g = the  acceleration  of  gravity 

(2)  The  rotational  inertia  of  the  concentric  solid 
cylindrical  components  of  the  passive  end  platen  and  its 
attachments  is  given  by 
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1 


u 


where:  w ^ = weight  of  the  i-th  solid  cylindrical  component 

d ,=  diameter  of  the  i-th  solid  cylindrical  component 
n = number  of  solid  cylindrical  components 


(3)  Transducers  and  other  masses  attached  to  this  platen 
can  be  accounted  for  by 


where:  vr  = weight  of  the  i-th  component 

rA = distance  from  the  platen  axis  to  thQ  center  of  mass 
for  the  i-th  component 

n = number  of  components  attached  to  passive  end  platen 
and  not  covered  in  the  determination  of  (Jp>i 


(4)  The  total  rotational  inertia  tor  the  passive  end  is 


given  by 


JP  = 1JP>1  * <JP>2 


c.  Active  End  Platen  Hass  and  Rotational  Inertia.  - (1)  The 
mass  , H A,  and  rotational  inertia,  JA , of  the  active  end  platen 
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shall  be  determined  with  all  transducers  and  rigid  attachments, 
including  attached  portions  of  the  vibration  excitation  device, 
securely  in  place. 

(2)  If  all  components  included  above  have  simple  geometry, 
then  the  equations  of  5 b.  may  be  used  to  obtain  the  mass  and 
rotational  inertia. 

(3)  An  alternate  procedure  involves  the  use  of  a metal 
calibration  rod  of  known  torsional  and  axial  stiffness.  One  end 
of  the  rod  shall  be  rigidly  fixed  and  the  other  end  shall  be 
rigidly  fastened  to  the  active  end  platen.  Since  it  may  be  very 
difficult  to  fasten  the  calibration  rod  to  the  platen  without 
adding  mass  or  rotational  inertia,  it  is  recommended  that  the 
calibration  rod  be  permanently  fastened  by  welding,  etc.,  to  an 
auxiliary  platen  identical  in  all  respects  to  the  one  to  be  used 
in  testing.  The  axial  and  torsional  stiffness  of  the  calibration 
rod  should  be  chosen  such  that  the  system  resonant  frequency  with 
the  calibration  rod  in  place  is  near  the  middle  of  the  range  of 
system  resonant  frequencies  anticipated  for  soil  testing. 
Several  calibration  rods  may  be  necessary  to  account  for 
different  specimen  sizes. 

(4)  Hith  the  calibration  rod  in  place,  determine  the  low 
amplitude  system  resonant  frequencies  for  longitudinal  vibration, 
(f  ro<j  and  for  torsional  vibration,  (f  rotj  )T  • The  mass  of  the 
active  end  platen  system  is  calculated  from 


"A  " 


<Krod>L 


(2^)  *2  [ (frod)L*2  - f oL  • 2 ] 
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calibration  rod 


where:  (Krod)L  = longitudinal  stiffness  of 
[ <AE)/L] 


A = calibration  rod  cross-sectional  area 
E = young's  modulus  for  calibration  rod  material 
L = length  of  calibration  rod 
foL  = apparatus  longitudinal  resonant  frequency  as 
described  in  5 d. 


The  rotational  inertia  of  the  active  end  platen  system  is 
calculated  from 


JA  = 


_<Wt. 


( 2ir)  4 2 [ (frod)x*2  - foX42] 


where: 


<Krod>X 


torsional  stiffness  of  calibration  rod  [ (IpG)/L] 
polar  moment  of  inertia  of  the  calibration  rod 


[ (it  d*4)/32  ] 


d = calibration  rod  diameter 

G = shear  modulus  for  calibration  rod  material 
foX  = apparatus  torsional  resonant  frequency  as 
described  in  5 d. 


The  above  equations  assume  that  the  mass  and  rotational  inertia 
of  the  calibration  rods  are  much  less  than  the  corresponding 
values  for  the  active  end  platen  system. 

(5)  A second  alternate  procedure  is  to  couple  the  metal 
calibration  rod  to  the  platens  in  place  of  the  specimen  and  then 
use  the  procedures  in  Section  8 to  backfigure  the  active  end 
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inertias  from  the  known  moduli  of  the  rod. 

d.  Apparatus  Resonant  Frequencies  and  Spring  Constants.  - (1) 
Apparatus  resonant  frequencies  and  spring  constants  are  defined 
only  for  those  apparatus  that  have  springs  attached  to  the  active 
end  platen  system  (See  7d.) . To  determine  the  resonant 
frequencies  , set  up  the  apparatus  complete  with  active  end 
platen  and  o-rings  but  no  specimen.  Vibrate  at  low  amplitude  and 
adjust  the  frequency  of  vibration  until  the  input  force  is  in 
phase  with  the  velocity  of  the  active  end  platen  system.  For 
longitudinal  vibration,  this  apparatus  resonant  frequency  is  f oL 
and  for  torsional  vibration  it  is  f oT  . 

(2)  The  longitudinal  and  torsional  apparatus  spring 

constants  (KgL  , KgT),  may  be  calculated  from 

KSL  * »A 

kST  * l^£or)(2  JA 

where  HA  and  JA  are  defined  in  5 c. 

e.  Apparatus  Damping  Factors.  - (1)  To  measure  the  damping 

constants  for  the  apparatus,  attach  the  same  masses  as  used  for 
the  determination  of  apparatus  resonant  frequencies.  For 

apparatus  without  springs  attached  to  the  active  end  platen, 
insert  the  calibration  rod  described  in  5 c.  With  the  apparatus 
vibrating  at  the  resonant  frequency,  cut  off  the  power  to  the 
excitation  device  and  record  the  decay  curve  for  the  vibration  of 


6 T = logarithmic  decrement  for  torsional  motion 

f.  Force/Current  and  Torque/Current.  - (1)  For  apparatus 
without  springs  attached  to  the  active  end  platen,  insert  the 
calibration  rod  as  described  in  5 c.  Determine  the  resonant 
frequency  of  this  single  degree-of-freedom  system  consisting  of 
the  active  end  platen  and  apparatus  spring  or  calibration  rod  by 
use  of  the  same  procedure  as  described  in  7 d. 

(2)  Then  set  the  frequency  to  0.707  times  the  resonant 
frequency  and  apply  sufficient  current  to  the  vibration 
excitation  device  so  that  the  vibration  transducer  output  to  the 
readout  device  has  a signa 1-to-noi se  ratio  of  at  least  10.  Read 
and  record  the  output  of  both  active  end  vibration  transducer  and 
current  measuring  instrument. 

(3)  Next  set  the  frequency  to  1.414  times  the  system 
resonant  frequency  and  obtain  the  vibration  transducer  and 
current  instrument  readings  in  a similar  fashion  to  those  at 
0.707  times  the  resonant  frequency. 

(4)  Calculate  and  C 2 from 

( VTCF)  ( TO  1 ) 

C1  4 

2 (CR  1) 

(VTCF)  , (T02) 

^ A 


CR2 

where:  (VTCF)A  = Active 

end  vibration 

transducer 

ca li brat  ion 

factor 

( LCF  a or  RCFa) 

depending 

on  whether 
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vibration  is  longitudinal  or  torsional. 

T0 1 = Active  end  transducer  output  at  0.707  x resonant 
frequency 

CR 1 = Current  instrument  reading  at  0.707  x resonant 
frequency 

T02  = Active  end  transducer  output  at  1.414  x resonant 
f requency 

CR2  = Current  instrument  reading  at  1.414  x resonant 
f requency 

C ^ and  c 2 should  agree  within  10  percent. 

(5)  By  use  of  and  C2  from  longitudinal  vibration,  the 
Force/Current  calibration  factor,  FCF , is  obtained  from 

FCF  = 0.  SfCj^  ♦ C2)  K SL 

where:  Kg^  = apparatus  spring  constant  (or  for  apparatus  without 
springs,  the  calibrating  rod  spring  constant)  for 
longitudinal  motion 

(6)  By  use  of  C^  and  C 2 from  torsional  vibration,  the 
Torque/Current  calibration  factor,  TCF,  is  obtained  from 

TCF  = 0. 5 (Cj  + C2)  K ST 


where:  K ST  = apparatus 

spr 

ing  constant 

(or 

for  ap 

paratus  without 

springs. 

the 

calibrating 

rod 

spring 
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torsional 

B 

O 
r * 

ion. 
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TEST  .SPECIMENS 


a.  Central.  - These  methods  cover  only  the  special  specimen 
preparation  procedures  related  to  the  vibration  and  resonant 
column  technique.  Since  the  resonant  column  test  may  be 
conducted  in  conjunction  with  controlled  ambient  stresses,  the 
provisions  for  preparation  of  specimens  in  Appendix  X or  Appendix 
XI  may  be  applicable  or  may  be  used  as  a quide  in  connection  with 
other  methods  of  application  and  control  of  ambient  stresses. 

b.  Specimen  Size.  - (1)  Specimens  shall  be  of  uniform 
circular  cross  section  with  ends  perpendicular  to  the  axis  of  the 
specimen . 

(2)  Specimens  shall  have  a minimum  diameter  of  33  am  (1.3 
in.)  and  the  largest  particle  contained  within  the  test  specimen 
shall  be  smaller  than  one-tenth  of  the  specimen  diameter  except 
that  for  specimens  having  a diameter  of  70  mm  (2.«  in.)  or 
larger,  the  largest  particle  size  shall  be  smaller  than  one-sixth 
of  the  specimen  diameter.  If,  after  completion  of  a test,  it  is 
found  that  larger  particles  than  permitted  are  present,  indicate 
this  information  in  the  report  of  test  data  under  "Remarks." 

(3)  The  length-to-diameter  ratio  shall  be  not  less  than  2 
(this  may  be  changed  to  1 for  torsional  vibration  only)  nor  more 
than  7 f except  that  when  an  ambient  axial  stress  greater  than  the 
ambient  lateral  stress  is  applied  to  the  specimen,  this  ratio 
shall  be  between  2 and  3. 

(4)  Measure  the  length  and  diameter  of  the  specimen  using  a 
Vernier  caliper  or  other  suitable  device.  Determine  the  weight 
of  the  test  specimen.  For  determination  of  moisture  content 
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(Appendix  I),  secure  a representative  sample  of  the  cuttings  from 
undisturbed  specimens,  or  of  the  extra  soil  for  remolded 
specimens,  placing  the  sample  immediately  in  a covered  container. 

c.  End  Coupling  for  Torsion.  - (1)  For  torsional  motion, 
complete  coupling  of  the  ends  of  the  specimen  to  the  specimen  cap 
and  base  must  be  assured.  Complete  coupling  for  torsion  may  be 
assumed  if  the  mobilized  coefficient  of  friction  between  the  end 
platens  and  the  specimen  is  less  than  0.2  for  all  shear  strain 
amplitudes.  (2)  The  coefficient  of  friction  is 
approximately  given  by 

Coeff.  of  Frict.  = V G/^ 

where:  Y = shear  strain  amplitude  (see  Section  8m.) 

G = shear  modulus  (see  Section  81.) 

o ' = effective  axial  stress 
a 

When  this  criterion  is  not  met,  other  provisions  such  as  the  use 
of  razor  blade  vanes,  adhesives,  etc.  must  be  made  in  order  to 
assure  complete  coupling. 

(3)  The  effectiveness  of  the  coupling  provisions  shall  be 
evaluated  by  testing  two  specimens  of  the  same  material  but  of 
different  lengths.  The  lengths  of  these  specimens  shall  differ 
by  at  least  a factor  of  1.5.  The  provisions  for  end  coupling  may 
be  considered  satisfactory  if  the  values  of  the  shear  modulus  for 
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these  two  specimens  of  different  length  do  not  differ  by  lore 
than  10  percent. 

d.  Specimen  Stiffness.  - (1)  The  upper  limit  of  specimen 
stiffness  will  depends  on  the  specific  apparatus  being  used. 
Refer  to  the  recommendations  of  the  manufacturer  or  to  Drnevich* 
for  methods  of  estimating  limits  of  specimen  stiffness. 

e.  Use  of  Side  Drains.  - (1)  For  tests  on  saturated  or  nearly 
saturated  soils  where  the  specimen  will  be  consolidated,  filter 
paper  strips  (Whatman  No.  54  or  similar)  may  be  used  around  the 
perimeter  of  the  specimen  to  speed  up  the  consolidation  and  pore 
pressure  equalization  processes.  Except  for  very  soft  specimens 
tested  at  low  effective  confining  stresses,  these  strips  should 
have  negligible  effects  on  the  test  results. 

(2)  Since  the  effect  of  side  drains  is  inversely  related 
to  specimen  diameter,  the  testing  of  larger  diameter  specimens 
can  be  used  to  check  on  the  effects  of  side  drains. 

(3)  For  side  drains  to  be  effective,  they  should  extend 
over  the  porous  stones  or  discs  at  each  end  of  the  specimen. 


7.  PROCEDURE 

a.  General.  - (1)  Resonant  column  testing  is  quite  similar  to 
triaxial  testing  except  that  vibrations  at  the  system  resonant 
frequency  are  applied  to  the  specimen  instead  of  the  deviator 
stress  that  causes  failure.  Consequently,  the  procedure  for 


* Drnevich,  V. P. , "Resonant  Column  Testing  - Problems  and 
Solutions,"  Proc.  of  the  ASTH  Symposium  on  Dynamic  Soil  and  Rock 
Testing  in  the  Field  and  Laboratory  for  Seismic  Studies,  STP  654, 
ASTH,  June,  1977. 
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preparation  of  specimens  in  Section  4.  of  Appendix  X.  is 
applicable  and  the  data  sheet  in  Plate  X-1  can  be  used  to  record 
the  specimen  data.  Because  the  application  of  vibrations  is 
relatively  nondestructive,  each  sample  can  be  tested  at  two  or 
three  confining  stresses  (staged  testing)  , and  thus  triplicate 
samples  are  not  usually  needed. 

(2)  Once  the  specimen  is  placed  in  the  apparatus,  a 
confining  stress  is  applied.  Depending  on  the  insitu  conditions 
being  simulated,  the  specimen  may  be  dry,  partially  saturated,  or 
saturated  by  means  of  a back  pressure.  Prior  to  the  application 
of  large  amplitude  vibration,  the  specimen  may  be  undrained  or 
consolidated.  (Low  amplitude  vibrations  where  strain  amplitudes 
are  on  the  order  of  0.001  percent  or  less  are  completely 
nondestructive  and  can  be  applied  at  any  time.) 

(3)  After  the  confining  stress  has  been  applied  (including 
completion  of  saturation  and  consolidation*  if  they  are  chosen) , 
small  amplitude  vibrations  at  the  system  resonant  frequency  are 
applied  and  data  are  taken.  Next,  the  amplitude  of  vibration  is 
increased  by  30  percent  to  50  percent,  the  frequency  is  adjusted 


* It  has  been  shown  by  D.  G.  Anderson  and  R.  D.  Woods 

( "Comparison  of  Field  and  Laboratory  Shear  Moduli,"  Pro- 
ceedings of  the  Geotechnical  Engineering  Specialty  Con- 
ference on  Insitu  Measurement  of  Soil  Properties,  ASCE, 
Raleigh,  N.C. , Vol.  I,  June  1975,  pp  69-92)  that  for  fine 
grained  soils,  the  amount  of  time  in  secondary  compression 
can  significantly  affect  test  results.  To  obtain  a measure 
of  this  effect,  it  is  recommended  that  low  amplitude  vibra- 
tion be  applied,  the  frequency  be  adjusted  to  resonance, 
and  that  sets  of  data  be  taken  at  convenient  intervals  of 
time  during  the  consolidation  process  and  to  at  least  one 
log  cycle  of  time  into  secondary  compression.  After  this, 
testing  at  larger  strain  amplitudes  can  commence. 
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to  resonance  and  another  set  of  data  are  taken.  This  process  is 
repeated  until  the  desired  amplitude  is  achieved  or  the  vibration 
capacity  of  the  apparatus  is  reached.  Then,  the  amplitude  of 
vibration  is  reduced  to  the  level  used  at  the  beginning  of  the 
test,  the  frequency  is  adjusted  to  resonance,  and  another  set  of 
data  is  taken. 

(4)  The  amount  of  disturbance  to  the  specimen  caused  by  the 
high  amplitude  vibration  can  be  established  by  comparing  the  low 
amplitude  data  obtained  before  and  after  high  amplitude 
vibration.  For  fine  grained  soils,  even  when  tested  with 
drainage  allowed,  excess  pore  pressures  may  be  created  by  the 
high  amplitude  vibration.  In  this  situation,  it  is  adviseable  to 
allow  the  excess  pore  pressures  to  dissipate  (specimen 
consolidates  or  swells).  Additional  sets  of  data  at  low 
amplitude  vibration  may  he  taken  during  this  process.  If  it  is 
apparent  from  the  data  that  only  a small  amount  of  disturbance 
was  caused  by  the  high  amplitude  vibration,  another  confining 
stress  may  be  applied  to  the  same  specimen  and  the  process 
outlined  above  can  be  repeated  for  the  new  confining  stress. 

b.  Procedure  for  Placement  of  Specimen  in  Apparatus.  - (1) 
The  exact  procedure  to  be  followed  during  test  setup  will  depend 
on  the  apparatus  and  electronic  equipment  used  and  on  the  methods 
used  for  application,  measurement,  and  control  of  the  ambient 
stresses.  However,  the  specimen  shall  be  placed  in  the  apparatus 
by  procedures  that  will  minimize  the  disturbance  of  the  specimen. 

(2)  Particular  care  must  be  exercised  when  attaching  the 
end  platens  to  the  specimen,  when  applying  filter  paper  strips. 
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membranes,  and  O-rings,  and  when  attaching  the  vibration 
excitation  device  to  the  platens.  A temporary  support  device  for 
the  platens  as  discussed  in  Section  4f  nay  be  needed. 

(3)  For  the  sake  of  good  coupling  between  the  specimen  and 
the  platens,  the  use  of  filter  paper  discs  at  the  ends  of  the 
specimen  is  not  recommended.  The  filter  paper  strips  on  the 
perimeter  of  the  specimen,  if  used,  should  extend  beyond  the  ends 
of  the  specimen  and  over  the  porous  discs  attached  to  the 
platens. 

(4)  When  it  is  desired  to  saturate  cohesive  soil  specimens, 
the  saturation  process  is  greatly  facilited  if  a very  small 
diameter  flexible  tube  is  inserted  between  the  top  platen  and  the 
membrane  before  the  top  0-rings  are  placed.  A seal  is  formed  by 
use  of  a liberal  amount  of  vacuum  grease.  A vacuum  is  then 
introduced  to  the  space  between  the  specimen  and  membrane  by 
means  of  this  tube.  Finally,  the  tube  is  gently  withdrawn  and  the 
top  O-rings  applied. 

(5)  For  cases  where  ambient  isotropic  stresses  are  to  be 
applied  to  the  membrane  enclosed  specimen,  liquid  or  air 
confining  media  may  be  used  for  dry  or  partially  saturated 
specimens.  For  tests  where  complete  saturation  is  important,  a 
liquid  confining  media  such  as  deaired  water  or  mineral  oil 
should  be  used.  where  the  vibration  excitation  device  is  located 
within  the  pressure  chamber,  an  air-liquid  interface  in  the 
chamber  just  below  the  device  is  acceptable  as  long  as  the  liquid 
covers  the  entire  membrane  that  encloses  the  specimen. 

(6)  The  system  for  providing  confining  stress  and  back 


pressures  to  the  speciaen  is  identical  to  that  given  in  Fig.  16 
of  Appendix  X.  except  that  a pressure  Line  must  be  connected 
between  valve  C on  the  chamber  pressure  reservoir  and  valve  B on 
the  apparatus  chamber  lid.  During  the  process  of  filling  the 
chamber  to  the  desired  level  with  confining  liquid,  valve  A is 
open.  After  this,  valve  A is  closed  and  valves  B and  C remain 
open . 

(7)  If  undisturbed  specimens  are  being  tested  and 
saturation  by  means  of  back  pressure  is  desired,  it  is 
recommended  that  consolidaton  be  done  before  back  pressuring. 
This  will  provide  a less  disturbing  effective  stress  history  to 
the  specimen. 

c.  Electric  Equipment  Connections.  - (1)  Connect  the 
vibration  excitation  device  to  the  sine  wave  generator  (with 
amplifier,  if  required) . The  power  supplied  to  the  vibration 
excitation  device  should  be  very  low  in  order  not  to  exceed  the 
amplitude  of  vibration  prescribed  in  7a. 

(2)  Connect  the  vibration  transducers  to  the  readout 
instruments  for  the  type  of  motion  (longitudinal  or  torsional)  to 
be  applied.  Adjust  the  readout  instruments  according  to  the 
instruction  manuals  for  these  instruments. 

d.  Measurement  of  Resonant  Frequency.  - (1)  The  procedure  for 
measuring  system  resonant  frequency  is  the  same  for  both 
longitudinal  and  torsional  vibration  except  that  the 
longitudinal  motion  transducer  is  used  tor  longitudinal  motion 
and  the  rotational  motion  transducer  is  used  for  torsional 
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motion 


(2)  If  the  passive  end  is  fixed  or  if  P>100  (See  Section  Be 
for  definition  ot  P) , motion  of  the  active  end  platen  is  used  to 
establish  resonance.  Otherwise,  motion  of  the  passive  end  platen 
is  used. 

(3)  tiith  the  power  as  low  as  practical,  increase  the 
frequency  of  excitation  from  a very  low  value  (e.g.  10  Hz.)  until 
the  resonant  frequency  is  obtained.  The  phase  relationship 
describing  resonance  can  be  established  by  observing  the 
Lissajous  figure  formed  on  an  x-y  oscilloscope  with  the  voltage 
proportional  to  the  driving  current  applied  to  the  horizontal 
amplifier  and  the  output  from  the  transducer  applied  to  the 
vertical  amplifier.  If  a velocity  transducer  is  used  for 
vibration  measurement,  the  system  resonant  frejuency  occurs  when 
the  figure  formed  is  a straight,  sloping  line.  If  a displacement 
or  acceleration  transducer  is  used,  the  frequency  should  be 
adjusted  to  produce  an  ellipse  with  axes  vertical  and  horizontal. 

(4)  For  the  case  where  the  motion  of  the  passive  end  platen 
is  being  used  to  establish  resonance,  the  resonant  frequency  is 
the  second  lowest  frequency  for  which  the  Lissajous  figure  either 
forms  a straight  sloping  line  (velocity  transducer  being  used)  or 
an  ellipse  with  axes  vertical  and  horizontal  (accelerometer  or 
displacement  transducer  being  used) . When  a velocity  transducer 
is  used  for  this  case,  the  slope  of  the  line  has  the  opposite 
sign  for  the  second  lowest  frequency  as  it  does  for  the  lowest 
frequency. 

(5)  It  is  recommended  that  the  frequency  be  measured  with  a 
digital  electronic  frequency  meter  and  be  recorded  to  at  least 
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three  significant  tigures.  The  system  resonant  frequency  for 
longitudinal  motion  shall  be  designated  f and  that  for  torsional 
motion  shall  be  designated  fT . 

e.  Measurement  of  Strain  Amplitude.  - (1)  The  strain 
amplitude  measurements  shall  only  be  made  at  the  system  resonant 
frequencies.  Thus,  for  a given  current  applied  to  the  excitation 
device,  the  vibratory  motion  transducer  outputs  recorded  at  the 
system  resonant  frequency  give  sufficient  information  to 
calculate  strain  amplitude.  To  increase  or  decrease  strain 
amplitude,  the  current  to  the  vibration  excitation  device  must  be 
increased  or  decreased.  After  making  a change  in  current 
applied  to  the  vibration  excitation  device,  the  procedure  of  7 d 
must  be  followed  to  establish  the  corresponding  system  resonant 
frequency  before  the  transducer  outputs  can  be  used  to  establish 
the  new  strain  amplitude  value. 

f.  Measurement  of  System  Damping,  - (1)  Associated  with  each 
strain  amplitude  and  system  resonant  frequency  is  a value  of 
damping.  Two  methods  are  available  for  measuring  system  damping: 
the  steady-state  vibration  method  and  the  amplitude  decay  method. 
The  procedures  for  both  methods  are  independent  of  whether 
longitudinal  or  torsional  motion  is  under  consideration. 

(2)  For  the  steady  state  method,  the  active  end  or  the 
passive  end  vibration  transducer  output  (depending  on  which  end 
is  used  to  establish  resonance)  and  the  current  applied  to  the 
vibration  excitation  device  must  be  measured  at  each  resonant 
frequency.  The  calculations  are  outlined  in  8 k. 

(3)  For  the  free  vibration  method,  with  the  system 
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vibrating  at  the  system  resonant  frequency,  cut  oft  the  pover  to 
the  vibration  excitation  device  and  record  the  transducer  output 
of  the  transducer  used  in  establishing  resonance  as  a function  of 
time.  This  gives  tne  decay  curve  for  free  vibration. 

g.  Explanation  of  Data  Sheet.  - (1)  Data  acquired  during  the 

performance  of  a test  is  recorded  directly  on  the  data  sheet 
given  in  Plate  1,  Columns  (2)  through  (7)  of  this  sheet  can 

be  used  to  record  consolidation  and/or  saturation  data.  Low 
amplitude,  nondestructive  vibration  data  may  also  be  recorded 
during  the  consolidation  and  saturation  processes. 

(2)  Data  in  column  (2)  is  the  elapsed  time  from  the 
beginning  of  the  test  or  from  the  initiation  of  consolidation. 

(3)  Data  in  columns  (3)  and  (4)  are  the  readings  from  the 
chamber  pressure  gage  and  the  pore  pressure  transducer, 
respectively  (See  Fig.  16,  Appendix  X).  They  may  be  in 
engineering  units  or  they  may  be  readings  to  which  calibration 
factors  must  be  applied.  For  the  latter  case,  the  calibration 
factors  for  each  should  be  noted  in  the  "Remarks"  section  at  the 
bottom  of  the  sheet.  If  a differential  pressure  transducer  (one 
side  measuring  cell  pressure  and  the  other  pore  pressure)  is 
used,  record  the  cell  pressure  from  the  pressure  gage  in  column 

(3)  and  the  differential  pressure  in  column  (4) . Relabel  column 

(4)  to  note  the  use  of  the  differential  pressure  transducer. 

(4)  Column  (5)  is  used  to  record  the  static  axial  load  that 
is  applied  in  addition  to  the  ambient  confining  stress. 

(5)  Columns  (6)  and  (7)  are  for  specimen  length  change 
measurement  readings  and  pore  volume  change  measurement  readings. 
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(6)  Place  an  "L"  or  a "T"  in  column  (8)  depending  on 
whether  longitudinal  or  torsional  excitation  is  being  applied  to 
the  specimen. 

(?)  Data  in  columns  (9)  through  (12)  are  associated  with 
the  vibration  of  the  specimen.  Column  (9)  is  used  for  the 
reading  associated  with  the  amount  of  force  or  torque  applied  to 
the  excitation  device  as  discussed  in  Sections  4c,  5e,  and  7f. 
The  outputs  of  the  active  end  and  passive  end  vibration 
transducers  are  recorded  in  columns  (10)  and  (11),  respectively. 
Finally,  the  system  resonant  frequency  (See  paragraph  d.  above) 
is  recorded  in  column  (12)  . 

(8)  If  the  free  vibration  method  is  also  used  to  measure 
damping,  each  record  of  amplitude  decay  with  time  must  be  marked 
to  associate  it  with  its  corresponding  line  in  the  data  sheet. 

8.  CALCULATIONS 

a.  General.  - (1)  Calculations  require  the  apparatus 
calibration  factors  (Section  5.)  and  the  physical  dimensions  and 
weight  of  the  specimen.  In  addition,  for  each  ambient  stress, 
one  data  set  is  required  for  each  vibration  strain  amplitude.  A 
data  set  consists  of  the  information  on  one  line  of  the  data 
sheet  given  in  Plate  1.  Additionally,  it  may  include  the 
free  vibration  amplitude  decay  curve. 

(2)  The  data  reduction  process  requires  several  steps,  the 
first  of  which  is  to  determine  the  specimen  volume,  length, 
diameter,  and  weight  for  each  data  set  (line  on  the  data  sheet) . 
Record  the  proper  values  of  these  specimen  characteristics  in 
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columns  (1)  through  (5)  of  the  Intermediate  Calculations  sheet 
given  in  Plate  2. 

(3)  Items  in  columns  (6)  through  (It)  of  this  sheet  are 
calculated  according  to  paragraphs  b.  through  g.  below.  Items  in 
columns  (8)  through  (11)  are  input  parameters  for  the  computer 
program  given  in  Section  1 1. 

(4)  If  the  computer  program  is  used,  values  in  columns  (12) 
and  (13)  are  calculated  by  the  program  and  the  information  in 
column  (14)  is  not  needed.  Alternatively,  items  in  columns  (12) 
through  (14)  may  be  determined  for  most  cases  from  equations  and 
graphs  discussed  in  paragraphs  h.  through  j. 

(5)  Information  for  the  Pinal  Calculations  Sheet  (Plate 

3) , columns  (2)  through  (5)  is  available  from  the  Data  Sheet 

(Plate  1).  Column  (2)  is  the  mean  effective  confining  stress 

(average  of  the  three  principal  effective  stresses).  Column  (3) 
is  the  ratio  of  axial  effective  stress  to  lateral  effective 
stress.  Column  (4)  is  simply  carried  over  from  column  (5)  of  the 
Data  Sheet.  Columns  (5)  through  (9)  are  determined  according  to 
procedures  given  in  paragraphs  k.  through  o.  below. 

f>.  Soil  Hass  Density.  - (1)  The  soil  mass  density,  p , is 
given  by 

v 

p = — 

Vg 

where?  H = total  weight  of  the  specimen  (Col.  (5) , Intermediate 
Calculations  Sheet) 

V = volume  of  the  specimen  (Col.  (2),  Intermediate 
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Calculations  sheet) 
g = acceleration  of  gravity 


Record  the  soil  mass  density  in  column  (6)  of  the  Intermediate 
Calculations  Sheet. 

c.  Specimen  Rotational  Inertia.  - (1)  The  specimen  rotational 
inertia  about  the  axis  of  rotation  is  given  by 

Wd|2 

8g 

where:  d = diameter  of  the  specimen  (Col.  (4),  Intermediate 

Calculations  Sheet) 


Record  the  value  of  specimen  rotational  inertia  in  column  (7)  of 
the  Intermediate  Calculations  Sheet. 

d.  Active  End  Inertia  Factors.  - (1)  The  active  end  inertia 
factor  for  longitudinal  motion,  T , is  calculated  from 


where:  = mass  or  the  active  end  platen  system  as  calculated 

in  5c. 

f oL  = apparatus  resonant  frequency  for  longitudinal  motion 
from  Sd.  (For  apparatus  without  springs  attached  to 
= the  active  end  platen,  this  term  is  zero.) 


f 


L 


= system  resonant  frequency  for  longitudinal  motion  (Col. 
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( 12)  , Data  Sheet) 


(2)  The  active  end  inertia  factor  for  torsional  motion,  TT  , 
is  given  by 


where:  JA  = rotational  inertia  of  the  active  end  platen  system 

as  calculated  in  5c. 

J = specimen  rotational  inertia  as  calculated  in  8c. 

= apparatus  resonant  frequency  for  torsional  motion 
from  5d.  (For  apparatus  without  springs  attached  to 
the  active  end  platen,  this  factor  is  zero.) 
fT=  system  resonant  frequency  for  torsional  motion  (Col. 
( 12) , Data  Sheet) 

Record  the  values  of  or  TT  in  column  (8)  of  the  Intermediate 

Calculations  Sheet. 

e.  Passive  End  Inertia  Ratios.  - (1)  For  longitudinal  motion, 
the  passive  end  inertia  ratio,  , is  given  by 


where:  M p = mass  of  the  passive  end  platen  system  as  described 

in  5b. 
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(2)  For  torsional  motion,  the  passive  end  inertia  ratio,  PT 
, is  given  by 


J 


where:  Jp  = rotational  inertia  of  passive  end  platen  system  as 

calculated  in  5b. 

For  the  special  case  where  the  passive  end  of  the  specimen  is 
rigidly  fixed,  PL  and  PT  are  equal  to  infinity.  For  most 
situations,  the  passive  end  may  be  assumed  to  be  rigidly  fixed  if 
the  value  of  P is  greater  than  10,000  . Record  passive  end 
inertia  ratio  in  column  (9)  of  Intermediate  Calculations  Sheet. 

f.  Apparatus  Damping  Factors.  - (1)  For  longitudinal  motion, 

the  apparatus  damping  factor,  ADFL , is  calculated  from 

ADFl  = ADCl/[  2TrfL(W/g)  ] 

where:  ADCL=  apparatus  damping  coefficient  for  longitudinal 
motion  as  described  in  5e. 

(2)  For  torsional  motion,  the  apparatus  damping  factor,  ADFT 
, is  calculated  from 

ADF  = ADC  /[  2irf  J ] 

T T T 

where:  ADC^  apparatus  damping  coefficient  for  torsional  motion 


\ 
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as  described  in  5e 


The  apparatus  damping  factors  are  recorded  in  column  (10)  of  the 
Intermediate  Calculations  Sheet. 

g.  Magnification  Factors.  - (1)  These  factors  are  used  in 

calculating  damping.  For  longitudinal  motion,  the  magnification 
factor  is  calculated  from 

MMF  L = [ (LCF)  (LTO)  / (FCF)  (C^)  ] (W/g)  (2frfL)  *2 

where:  LCF  = longitudinal  motion  transducer  calibration  factor 
for  the  transducer  used  in  establishing  resonance 
LTO  = longitudinal  motion  transducer  output  (Col.  (10)  or 
(11),  Data  Sheet  depending  on  which  end  is  used  to 
establish  resonance.) 

FCF  = force/current  factor  from  5f. 

CR  = current  reading  to  longitudinal  excitation  system 
L 

(Col.  (9)  , Data  Sheet) 

(2)  For  torsional  motion,  the  magnification  factor  is 
calculated  from 

MMFt  = [ ( RCF)  (RTO)  / (TCF)  (CR  T)  ]J(2trfT)*2 

where:  RCF  = rotational  transducer  calibration  factor  for  the 
transducer  used  in  establishing  resonance 
RTO  = rotational  transducer  output  (Col.  (10)  or  (11), 
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Data  Sheet  depending  on  which  end  is  used  to 
establish  resonance.) 

TCF  = torque/current.  factor  from  5f. 

CR T = current  reading  to  torsional  excitation  system  (Col. 
(9) , Data  Sheet) 

Record  magnification  factors  in  column  (11)  of  Intermediate 
Calulations  Sheet. 

h.  Dimensionless  Frequency  Factor.  - (1)  The  dimensionless 
frequency  factor,  F,  is  used  in  calculating  modulus.  It  is  a 
function  of  system  factors  T,  P,  and  ADF  and  of  specimen  damping 
ratio,  D.  Values  of  F are  provided  by  the  computer  program  in 
Section  11. 

(2)  For  cases  where  ADF  is  zero  and  specimen  damping  ratio 
is  less  than  10%,  values  of  F can  be  obtained  from  Figure  2. 
Figure  2b  is  similar  to  Figure  2a  except  that  the  range  of  T is 
different.  This  figure  is  independant  of  which  end  of  the 
specimen  is  used  to  determine  resonance. 

(3)  For  the  cases  where  ADF  and  specimen  damping  ratio  are 
small  and  where  T and  P are  both  greater  than  10,  dimensionless 
frequency  factor  may  be  approximated  by 

F = V 1/T  ♦ 1 I/P 

This  simple  equation  is  especially  useful  for  checking  results 
obtained  from  Fig.  2 or  from  the  computer  program.  Record 
dimensionless  frequency  factors  in  column  (12)  of  the 
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Intermediate  Calculations  Sheet. 

i.  Strain  factor.  - (1)  The  strain  factor,  SF,  is  calculated 
by  the  program  in  Section  11.  or,  for  cases  of  ADF  = 0 and 
specimen  damping  equal  to  10  percent,  it  may  ba  obtained  from 
Fig.  3.  For  other  values  of  specimen  damping  ratio  and  ADF, 
values  of  from  Fig.  3.  are  only  approximately  correct.  (Note  that 
Fig.  3a  is  for  the  case  where  resonance  is  established  by  phase 
mesasurement  between  input  force  and  motion  at  the  active  end  and 
Fig.  3b  is  for  the  case  where  resonance  is  established  by  phase 
measurement  between  input  force  and  motion  at  the  passive  end.) 
Record  strain  factors  in  column  (13)  of  the  Intermediate 
Calculations  Sheet. 

j.  Amplification  Coefficient.  - (1)  When  manually  reducing 

data,  the  amplification  coefficient.  A,  is  needed  to  calculate 
damping.  This  coefficient  may  be  obtained  from  Fig.  4 for  cases 
where  apparatus  damping  factor,  ADF  (Col.  (10) , Intermediate 
Calculations  Sheet),  is  zero.  For  other  values  of  ADF,  the 
amplification  coefficient  is  only  approximate. 

(2)  The  following  equation  may  also  be  used  to  approximate 
the  amplification  coefficient  for  the  case  where  resonance  is 
determined  by  phase  measurement  between  input  force  or  torque  and 
motion  at  the  passive  end.  For  this  equation,  the  values  of  both 
T and  P must  be  greater  than  10. 


A = 2 (T  + P) 

Amplificaton  coefficient  is  recorded  in  column  (14)  of  the 
Intermediate  Calculations  Sheet. 

k.  Number  of  Vibration  Cycles.  - (1)  The  approximate  number  of 


J 


vibration  cycles,  N,  applied  for  each  data  set  may  be  determined 
by 


N = 60(tt  - t^)  f 

; 

where:  t = elapsed  time  in  minutes  for  data  set  "i" 

tj_j  = elapsed  time  in  minutes  for  data  set  "i-1" 

f = system  resonant  frequency  (Col.  (12) , Data  Sheet) 

I 

Both  tA  and  ti_1  are  from  column  (2)  of  the  Data  Sheet.  Number  of 
vibration  cycles  is  recorded  in  column  (5)  of  the  Final 
Calculations  Sheet. 

1.  Nodulr.  - (1)  The  rod  modulus  is  calculated  from 

E = P(2ttL)*2  (fL/FL)»2 

where:  P = specimen  mass  density  (Col.  (6),  Intermediate 

Calculations  Sheet) 

f^  = system  resonant  frequency  for  longitudinal  motion 
(Col.  (12),  Data  Sheet) 

P = dimensionless  frequency  factor  (Col.  (12), 
L 

Intermediate  Calculations  Sheet) 

(2)  The  shear  modulus  is  calculated  from 

G = P (2ftL)  *2  (f  /F  ) 1 2 
T T 


L 
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where:  = system  resonant  frequency  for  torsional  motion  (Col. 

(12),  Data  Sheet) 

Ft  = dimensionless  frequency  factor  (Col.  (12), 
Intermediate  Calculations  Sheet) 

Record  moduli  in  column  (6)  of  Final  Calculations  Sheet. 

a.  Strain  Amplitude.  - (1)  The  averaqe  rod  strain 

amplitude,  e , for  longitudinal  motion  shall  be  calculated  from 

e = (LCF) (LTO) (SF/L)  100* 


where:  LCF 


LTO 


SF 


L 


longitudinal  motion  transducer  calibration  factor 
for  the  transducer  used  in  establishing  resonance 
longitudinal  transducer  output  (Col. (10)  or  (11), 
Data  Sheet  depending  on  transducer  used  in 
establishing  resonance) 

strain  factor  (Col.  (13),  Intermediate  Calculations 
Sheet ) 

specimen  length  (Col.  (3),  Intermediate  Calculations 
Sheet ) 


(2) 


amplitude. 


For  torsional  motion,  the 
Y,  shall  be  calculated  from 


average 


shear  strain 


Y = (RCF)  (SF)  [ d/(3L)  ] 100% 


where:  RCF  = rotational  motion  transducer  calibration  factor  for 
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the  transducer  used  in  establishing  resonance 


RTO  = rotational  transducer  output  (Col.  (10)  or  (11), 
Data  Sheet  depending  on  transducer  used  in 
establishing  resonance) 

SF  = strain  factor  (Col.  (13),  Intermediate  Calculations 
Sheet) 

d = specimen  diameter  (Col.  (4) , Intermediate 

Calculations  Sheet) 

Strain  amplitudes  are  recorded  in  column  (7)  of  Final 
Calculations  Sheet. 

n.  Damping  Ratio  from  Steady  State  Vibration.  - (1)  If  the 
computer  program  in  Section  11  is  used,  the  damping  ratio  of  the 
specimen  is  established  as  part  of  the  output. 

(2)  Manual  calculation  of  damping  ratio  may  be  done  for 
cases  where  the  apparatus  damping  factor,  A DF , is  zero  or  may  be 
assumed  to  be  zero  (See  8f . for  definition  of  ADF) . 

The  damping  ratio,  D,  is  calculated  from 

D = 100*/[  A (MP1F)  ] 

where:  A = Amplification  Coefficient  (Col.  (14) , Intermediate 

Calculations  Sheet) 

MMF  = magnification  factor  (Col. 

Calculations  Sheet) 


4 1 


(11) 


Intermediate 


o.  Damping  Ratio  from  Free  Vibration.  - (1)  This  procedure  is 
theoretically  exact  for  apparatus  where  the  passive  end  can  be 
assumed  to  be  rigidly  fixed.  For  cases  where  the  passive  end  is 
not  rigidly  fixed,  irrespective  of  which  end  of  the  specimen  is 
used  in  establishing  resonance,  this  method  is  approximate.  The 
same  transducer  that  is  used  to  determine  resonance  must  be  used 
to  obtain  the  amplitude  decay  curve. 

(2)  For  the  case  where  resonance  is  established  by  use  of 
the  passive  end  transducer,  values  of  T and  P should  both  be 
greater  than  ten  when  amplitude  decay  is  used. 

(3)  For  apparatus  where  the  active  end  platen  is  restrained 
by  a spring,  a system  energy  ratio  must  be  calculated.  For  other 
apparatus,  this  factor  is  zero. 

(4)  For  longitudinal  motion,  the  system  energy  ratio  is 
calculated  from 


SL  = <V/W)<foLFL/fL>*2 


where:  F^  = dimensionless  frequency  factor  for  longitudinal 

motion  (Col.  (12),  Intermediate  Calculations  Sheet) 

and  tor  torsional  motion  from 


= (J  /J) 

A 


<f 


oT  T 


'V 


♦ 2 


where:  FT  = dimensionless  frequency  factor  for  torsional  motion 

(Col.  (12),  Intermediate  Calculations  Sheet) 
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Record  the  system  energy  ratios  in  column  (7)  of  the  Amplitude 
Decay  Data  Sheet. 


(5)  Compute  the  system  logarithmic  decrement  from  the  free 
vibration  (obtained  in  7f.)  from 

6S  = (1/n)  log  (A!/An+1) 

where:  = amplitude  of  vibration  for  first  cycle  after  power 

is  cut  off  (Col.  (3),  Amplitude  Decay  Data  Sheet) 

An+i  = amplitude  of  vibration  for  (n+1)th  cycle  of  free 
vibration  (Col.  (4) , Amplitude  Decay  Data  Sheet) 
n = number  of  free  vibration  cycles  which  must  be  10  or 
less.  (Col.  (5)  , Amplitude  Decay  Data  Sheet) 

(6)  Finally,  calculate  the  damping  ratio  from 

D = [ 6s(1  + S)  -S  6]100V(2ir) 

D L or  DT  , in  percent,  depending  on  whether 
vibration  is  longitudinal  or  torsional 
&sL  or  6gT  depending  on  whether  vibration  is 

longitudinal  or  torsional 

S^  or  ST  depending  on  whether  vibration  is 

longitudinal  or  torsional 

6^  or  6t  , the  apparatus  logarithmic  decrement  from 
5e. 
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where:  D = 


6 = 

s 


S = 


6 = 


Record  values  of  damping  ratio  in  both  column  (8)  of  the 
Amplitude  Decay  Data  sheet  and  in  column  (9)  of  the  Pinal 
Calculations  Sheet. 

p.  Accuracy  of  Calculations.  - (1)  The  items  listed  on  the 
Final  Calculations  Sheet  are  those  which  will  be  used  in  analysis 
or  design  and  consequently,  it  is  important  that  they  be 
accurate. 

(2)  In  the  foregoing  sections,  no  units  are  specified.  Any 
system  of  units,  whether  English,  metric,  or  SI,  will  work  as 
long  as  units  are  used  in  a consistent  fashion. 

(3)  All  numerical  items  in  the  Data  Sheet  (Plate  1) 
should  be  recorded  to  at  least  three  significant  figures.  In 
cases  where  apparatus  have  springs  attached  to  the  active  end 
platen,  and  where  the  system  resonant  frequencies  are  within  20 
percent  of  the  apparatus  resonant  frequency,  the  system  resonant 
frequency  should  be  recorded  to  at  least  four  significant 
figures.  A digital  frequency  meter  is  generally  required  for 
this  accuracy  and  when  the  system  resonant  frequencies  are  less 
than  100  Hz,  it  may  be  easier  and  more  accurate  to  measure  the 
resonant  period  and  invert  it  to  get  the  resonant  frequency. 

(4)  Items  in  columns  (2)  through  (7)  in  the  Intermediate 

Calculations  Sheet  (Plate  2)  should  be  calculated  and 

recorded  with  three  significant  figures.  Any  constants  or  units 
conversions  for  these  items  should  also  be  to  at  least  three 
significant  figures.  Items  in  columns  (8)  througn  (12)  are  all 
non-dimensional.  It  is  imperatative  that  the  parameters  used  in 
calculating  these  items  be  in  consistent  units.  Special 
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precaution  should  be  taken  when  determining  the  active  end 
inertia  ratio,  T,  in  column  (8)  and  subsequent  parameters  which 
are  functions  of  this  inertia  ratio,  when  the  system  resonant 
frequency  is  within  20  percent  of  the  apparatus  frequency.  For 
this  case,  it  is  also  adviseable  to  use  the  computer  program  to 
determine  the  dimensionless  frequency  factor,  F,  the  strain 
factor,  SF,  and  the  damping  ratio,  D. 

(5)  If  all  of  the  above  guidelines  are  observed,  the  data 
on  the  Final  Calculations  Sheet  (Plate  3)  should  be  reported 
to  three  significant  figures.  The  absolute  accuracy  of  these 
data  will  then  be  directly  related  to  the  absolute  accuracy  of 
the  calibrations. 

9.  REPORT 

a.  General.  - The  report  shall  include  characteristics  of  the 
apparatus,  specimen,  ambient  test  conditions,  and  the  results  for 
each  data  set. 

b.  Apparatus  Characteristics.  - The  following  apparatus 

characteristics  shall  be  included:  apparatus  name,  model  number 

and  serial  number;  active  end  and  passive  end  masses  and 

rotational  inertias  (HA,  Up,  JA  , J pj  ; longitudinal  and  torsional 

apparatus  resonant  frequencies  (foL  * foT^;  longitudinal  and 

torsional  apparatus  logarithmic  decrements  (5  ,5  ) ; the 

jfc  T 

force/current  and  torque/current  constants  (FCF,  TCF) ; and  the 
applicable  motion  transducer  calibration  factors  (LCFA  , LCFp  , RCF  A 
, RCF p ).  (Note  that  if  the  passive  end  is  fired,  inertias  and 
transducers  are  not  needed  for  the  passive  end.  Likewise,  if 
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only  one  type  of  notion,  longitudinal  or  torsional,  is  used,  then 
only  factors  and  inertias  for  that  type  need  be  given.) 

c.  Specimen  Characteristics.  - A visual  description  and 
origin  of  the  soil  shall  be  given  including:  name,  group  symbol, 
and  whether  undisturbed  or  remolded.  This  may  be  reported  on  the 
Specimen  Data  Sheet  (Plate  X-1). 

d.  Ambient  Stress  History.  - This  is  to  include  a description 
of  the  sequence  of  loading,  consolidating,  and  saturating  the 
speciaen  and  the  application  of  axial  loads.  For  staged  testing 
where  behavior  is  measured  at  more  than  one  effective  confining 
stress,  the  details  of  this  staging  should  be  described.  Plots 
of  specimen  volume  change  and/or  specimen  length  change  versus 
log  time  or  square-root-of-time  during  consolidation  should  also 
be  given. 

e.  Vibration  Test  Results.  - This  is  to  include  all  of  the 
applicable  information  on  the  data  and  calculation  sheets 
(Plates  1 through  -4)  for  each  data  set  including  those  associated 
with  consolidation  and  saturation.  For  each  stage  of  testing, 
plots  of  modulus  and  damping  (arithmetic  scales)  versus  vibratory 
strain  amplitude  (either  arithmetic  or  logarithmic  scale)  should 
also  be  given. 

10.  POSSIBLE  ERRORS 

a.  Specimen  Preparation.  - (1)  since  specimen  preparation  is 
identical  to  that  for  triaxial  testing,  the  same  sources  of  error 
discussed  in  Appendix  X.,  especially  in  Section  9 b.  are 
applicable  here. 
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b.  Apparatus  Assembly.  - (1)  The  most  common  source  of  error 
associated  with  apparatus  assembly  is  disturbing  the  specimen 
when  placing  the  end  platens  and  when  attaching  the  vibration 
excitation  device  to  the  platen.  A temporary  platen  support 
device  as  described  in  Section  4 g.  and  a temporary  support  for 
the  vibration  excitation  device  as  described  in  Section  4 f. 
should  be  used  to  avoid  disturbance. 

(2)  Complete  contact  between  the  ends  fo  the  specimen  is 
essential  for  coupling  in  torsion  and  for  accurate  rod  modulus 
and  rod  strain  values.  Careful  use  of  the  mitre  box  in  trimming 
the  specimen  ends  will  circumvent  the  problem  for  cohesive  soils 
and  careful  smoothing  and  leveling  of  the  specimen  top  will 
circumvent  this  problem  for  cohesionless  soils. 

(3)  When  stiff  to  very  stiff  soils  are  tested  and  when 
razor  blade  vanes  are  protruding  from  the  end  platens,  there  is  a 
possibility  that  the  blades  will  not  fully  penetrate  the 
specimen.  This  can  be  averted  in  most  cases  by  shortening  the 
distance  the  blades  protrude  from  the  porous  discs  and  by  keeping 
the  blades  razor  sharp.  For  extremely  stiff  soils,  the  blades 
should  not  be  used.  Capping  adhesives  serve  to  provide  the  proper 
contact  for  these  soils. 

(4)  Most  apparatus  rely  on  a coil-permanent  magnet  system 
for  vibration  excitation.  For  cases  where  the  active  end  platen 
does  not  have  springs  attached  to  it,  the  coils  may  touch  the 
magnets,  especially  as  the  specimen  swells  or  consolidates. 
Usually  there  are  provisions  for  adjusting  the  clearance  between 
the  coil(s)  and  magnet  (s)  without  having  to  remove  the  confining 
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stresses  or  otherwise  disturb  the  specimen.  Adjustments  should 
be  made  as  needed. 

(5)  Associated  with  the  vibration  excitation  device,  the 
transducers,  and  the  pore  water  control  are  various  cables  and 
tubes.  It  is  important  that  these  be  arranged  such  that  they  do 
not  affect  the  vibratory  characteristics  of  the  system. 

c.  Test  Procedure.  - (1)  The  definition  of  resonance  is 
different  for  different  apparatus  boundary  conditions. 
Particularly  where  the  motion  at  the  passive  end  is  used  to 
measure  resonance,  the  wrong  resonant  freguency  may  be  recorded 
as  the  system  resonant  frequency.  Also,  the  apparatus  may  have 
some  minor  spurious  resonances  near  the  system  resonant 
frequency.  These  should  be  known  beforehand  by  use  of  the 
calibrating  rod  and  should  not  be  mistaken  tor  the  system 
resonant  frequency. 

(2)  Another  source  of  error  in  the  apparatus  where  the 
passive  end  platen  is  not  rigidly  fixed,  is  the  use  of  the  wrong 
transducer  to  establish  resonance.  This  problem  can  be  averted 
by  clearly  labeling  all  leads  from  the  apparatus. 

(3)  At  all  amplitudes  of  vibration,  but  particularly  at 
higher  strain  amplitudes,  the  system  resonant  frequency  changes 
with  time  due  to  changes  occuring  within  the  specimen  from 
vibration  and/or  consolidation.  If  the  voltage  readings  from  the 
transducers  are  not  taken  while  the  system  is  at  resonance, 
significant  errors  in  strain  amplitude  and  damping  ratio  will 
occur. 

d.  Data  Seduction.  - (1)  Other  than  the  common  errors  of 
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calculation  and  graph  interpretation,  errors  due  to  units 
conversion  and  use  of  inconsistent  units  are  relatively  frequent. 
It  is  preferred  that  one  system  of  units  be  used  throughout. 
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(2)  Finally,  it  is  important  that  all  results  be  checked 
for  reasonableness.  Approximate  values  for  test  results  may  be 
obtained  by  use  of  the  equations  and  curves  by  Hardin  and 
Drnevich.*  Gross  deviations  from  expected  values  or  inconsistent 
results  for  a given  test  should  be  cause  for  rechecking  data  and 
calculations  and  possibly  redoing  the  test. 

11.  COMPUTE R PROGRAM  FOR  DATA  REDUCTION 

a.  General.  - (1)  All  resonant  column  data  associated  with 
describing  the  system,  the  amplitude  of  vibration,  and  resonant 
frequency  can  be  represented  by  five  nondimensional  terms  that 
are  independent  of  apparatus  type  and  whether  the  vibratory 
motion  is  longitudinal  or  torsional.  Four  of  the  terms  (T,  P, 
ADF,  and  MMF)  are  given  in  columns  (8)  through  (11)  of  the 
Intermediate  Calculations  Sheet.  The  remaining  term,  JPA, 
describes  which  end  (active  or  passive)  is  used  in  establishing 
resonance. 

(2)  The  FORTRAN  program  adapted  from  Drnevich  et  al.** 
calculates  three  nondimensional  parameters.  The  first  is  the 
dimensionless  frequency  factor,  F,  which  is  used  along  the  system 


* Hardin,  B.O.,  and  Drnevich,  V.P., "Shear  Modulus  and 
Soils:  Design  Equations  and  Curves,"  Journal  of 

Mechanics  and  Foundations  Division,  ASCE,  Vol.  98,  No. 
1972,  pp.  667-692. 

**Drenvich,  Hardin,  and  Shippy,  op.  cit. 


Damping  in 
the  Soil 
SM7,  July, 
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resonant  frequency  and  specimen  characteristics  to  calculate 
modulus.  The  second  is  the  strain  factor,  SF,  which  is  used 
along  with  vibration  amplitude  and  specimen  dimensions  to 
calculate  the  average  vibratory  strain  amplitude.  The  last 
parameter  is  the  damping  ratio  of  the  specimen,  in  percent,  which 
is  calculated  directly  by  the  program. 

(3)  The  program  utilizes  a Gaussian  elimination  procedure* 


to  solve 

four 

si multaneous 

equations. 

A double 

iteration 

procedure 

(one 

for  di mensionl ess  frequency 

factor  and 

one  for 

damping 

rat io) 

is  achieved  by 

using  two 

M ueller ' s 

iteration 

subroutines**. 

Each  iteration 

subrout ine 

requires 

an  error 

criterion 

(EPSD 

is  used  for  damping  ratio 

and  EPSF  is 

used  for 

dimensionless  frequency  factor). 

Both  may 

be  prescribed  Dy  the 

user  or  default  values  built  into  the  program  may  be  used. 

(4)  In  addition  to  the  output  described  above,  the  program 
outputs  a calculated  phase  angle  and  a calculated  magnification 
factor  that  may  be  used  for  the  purpose  of  checking  the  accuracy 
of  the  iterative  processes.  The  phase  angle  should  be  equal  to 
zero  and  the  calculated  magnification  factor  should  equal  the 
input  magnification  factor,  MHF,  which  is  based  on  measurements 
and  calibration  factors.  If  the  differences  between  the  measured 
and  calculated  values  of  flHF  are  more  than  two  percent  a decrease 
in  EPSD  should  be  made  until  the  difference  is  within  two 
percent.  Similarly,  if  the  phase  is  greater  than  0.001  radian,  a 


* System/360  Scientific  Subroutine  Package,  Version  III,  Fifth 
Edition,  18*,  1970. 

**  ibid. 
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decrease  in  EPSF  should  be  made.  Very  low  values  of  EPSD  and 
EPSF  should  be  avoided  because  they  cause  excessive  iteration 
cycles  and  coiputer  tine. 

b.  Input  Data  Format.  - (1)  The  first  card  sets  the 
dimensionless  frequency  factor  error  criterion,  EPSF,  the  maximum 
number  of  iterations  for  dimensionless  frequency  factor,  ITERF, 
the  damping  error  criterion,  EPSD,  and  the  maximum  number  of 
iterations  for  damping  ratio,  ITERD.  The  format  for  these  data 
are  2 ( FI 0. 0, 7X , 13) . The  default  values  of  EPSF  = 0.01,  ITERF  = 
40,  EPSD  = 0.0001,  and  ITERD  = 40.  To  use  the  default  values, 
put  a blank  card  in  for  the  first  card. 

(2)  Next,  there  is  one  card  for  each  data  set  (line  on 
Intermediate  Calculations  Sheet).  The  data  required:  T,  P,  ADF, 
anF,  and  J PA  comes  from  columns  (8)  through  (11)  of  the 
Intermediate  Calculations  Sheet  for  the  first  four  and  the  value 
of  J P A equals  0 when  the  active  end  is  used  for  establishing 
resonance  and  equals  1 when  the  passive  end  is  used.  The  format 
for  these  data  sets  is  4 FI  0.  0,  9X , 1 1.  Any  number  of  data  sets  may 
be  used  as  long  as  the  last  card  in  the  data  deck  is  a blank 
card . 

Sample  Input  Data  Set  for  Program 


1 2 3 A 5 

123  **56 78 90 123  *•*'678901  23 A6678901  23 A 667  P90 1 23 AS  67890  Data  Card  Col.  Nos. 


0.01 

050  0.0001 

050 

First  Card 

98.75 

A.  69 

.1053 

.6729 

1 

97.98 

A. 69 

.1070 

• 5A66 

l 

Data  Set  Cards 

96.30 

A. 69 

.1108 

.3985 

1 

9A.3A 

*••69 

.1150 

.3009 

1 

Last  Card  is  Blank 
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FILE:  RESCOL 


PROG 


A 


CONVERSATIONAL  MONITOR  SYSTEM 


. COMPUTER  PROGRAM  FOR  RESONANT  COLUMN  DATA  REDUCTION 


DEFINITIONS  OF  INPUT-OUTPUT  VARIABLES 

ADF  APPARATUS  DAMPING  FACTOR  ( ADF  > 0.) 

f)  SPECIMEN  DAMPING  RATIO  (0.01?  < D < 35?) 

EPSD  ERROR  CRITERION  FOR  D (DEFAULT  VALUE:  O.OOD1) 

EPSF  ERROR  CRITERION  FOR  F (DEFAULT  VALUE:  0.01) 

F FREQUENCY  FACTOR 

ITERD  MAXIMUM  NUMBER  OF  ITERATIONS  ALLOWED  FOR  D (DFFAULT 

VALUE:  AO) 

ITERF  MAXIMUM  NUMBER  OF  ITERATIONS  ALLOWEO  FOR  F (DEFAULT 

VALUE:  AO) 

JPA  INDICATOR  OF  END  WFfE  RE  MEASUREMENTS  WFRE  TAKEN: 

J PA  = 0 FOR  MEASUREMENTS  AT  THF  ACTIVE  END; 

JPA  = 1 FOR  MfcASUREMt  NTS  AT  THE  PASSIVE  END. 

MMF  MODIFIED  MAGNIFICATION  FACTOR  (MMF  > 0.) 

P PASSIVE-END  INERTIA  RATIO  (P  > O.I  IF  JPA  = Or  THEN 

P > 100.  AND  P > T)  • 

sf  strain  factor 

T ACTIVE-END  inertia  FACTOR  (T  > -10.) 


GIVEN  values  OF  T,  P,  ADF.  MMF,  AND  JPA,  THIS  PROGRAM  CALCU- 
LATES VALUES  OF  F,  D,  AND  SF  AND  WILL  PRINT  VALUES  OF  ALL  THESE 
PARAMETERS. 


RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 

RCP 


RCP 


DATA  INPUT  INSTRUCTIONS RCP 

RCP 

VALUES  OF  EPSF,  ITERF,  EPSD,  AND  HERD  MAV  BE  SPECIFIED  ON  THF  RCP 
FIRST  DATA  CARD  ACCORDING  TO  THE  FORMAT  2 ( F 10 .0, 7X , 13 ) . THE  SF  RCP 

»ARAMETtRS  ARE  REQUIRED  IN  THE  ITERATIVE  SOLUTION  PROCEDURE  TO  CON-  RCP 
1R0L  ACCURACY  AND  LIMIT  THE  NUMBER  OF  ITERATIONS.  THIS  DATA  CARD  RCP 
MAY  BE  LEFT  BLANK  IF  THE  USER  CHOOSES  TO  USE  THE  DEFAULT  VALUES  RCP 

SPECIFIED  BY  THE  PROGRAM.  (SEE  THE  LIST  OF  DEFINITIONS  OF  INPUT-  RCP 


OUTPUT  VARIABLFS.)  THE  DFFAULT  VALUES  WILL  GIVF  GOOD  R F SULTS  IN  MOSTRCP 
CASES  INVOLVING  SMALL  APPARATUS  DAMPING.  FOR  LARGE  APPARATUS  DAMPINGRCP 
UR  TU  CHECK  RESULTS  OBTAINED  WITH  THE  DEFAULT  VALUFS,  SMALLER  VALUES  RCP 
UF  ITERF  AND  ITERD  SHOULD  BE  USED.  (HOWEVER,  IN  SOME  CASES  OF  LARGE  RCP 


APPARATUS  DAMPING,  ACCURATE  CALCULATION  OF  0 IS  IMPOSSIBLE.)  RCP 

RCP 

EACH  SUBSEQUENT  OATA  CARD  SHOULD  CONTAIN  A VALUE  OF  EACH  OF  THE  RCP 
PARAMETERS  P,  T,  ADF,  MMF,  AND  JPA  ACCORDING  TO  THE  FORMAT  AF10.0,  RCP 
NX, II.  THUS,  EACH  OF  THFSE  CARDS  CORRESPONDS  TO  ONE  SET  OF  EXPERI-  RCP 
MENIAL  TEST  DATA.  SEE  THE  LIST  OF  DEFINITIONS  OF  INPUT-OUTPUT  RCP 

VAR1ABLLS  FOR  THF  RANGES  OF  THESE  VARIABLES  ALLOWED  BY  THE  PROGRAM.  RCP 
THE  DATA  CAROS  SHOULD  BE  TERMINATED  WITH  A BLANK  CARO.  RCP 

RCP 

RCP 

RCP 

DIMENSION  C ( A ) , Up ( 1 1 ) RCP 

REAL  MMF ,MMFCAL  RCP 


0001 

0002 

0003 

000  A 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 
0013 
001 A 

001  A 
0016 

0017 

0018 
001« 
0020 
0021 
0022 
0023 
002A 

0025 

0026 

0027 

0028 

0029 

0030 
00’ 1 

0032 

0033 
003A 

0035 

0036 

0037 

0038 

0039 
OOAO 
OOA  1 
00A2 
00A3 
OOAA 
OOAS 
00A6 
OOA  7 
00A8 
OOAO 
0050 
OOA  | 

0052 

0053 
005  A 
0055 
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FILE 


RESCOL 


PROG 


A 


CONVERSATIONAL  MONITOR  SYSTEM 


INTEGER  00  RCP 

COMMON  /VALS/P.T,ADF,MMF,JPA,AMP,0,F  /DLIM/OL.DR  RCP 

COMMON  /ALARM/KALARM  /E PS  I T/EPS F , I T ERF , EPSO , I TE RD  RCP 

COMMON  /PAR/PL,QL,C  RCP 

COMMON  /CRIT/MMFCAL, PHASE  RCP 

tXTEKNAL  DELAMP  RCP 

1000  FORMAT (4F10.0,9X,I1 > RCP 

1005  FORMAT  ( 2!  F10  .0  « 7X  » I 3)  ) RCP 

1110  FORMAT ( 49H  *****  WARNING  *****  POSSIBLY,  RCP 

1 23H  NOT  ENOUGH  ITERATIONS  ) RCP 

1115  FORMA T ( 50H  TO  OBTAIN  SPECIFIED  ACCURACY  FOR  D WITH  THE  ABOVE,  RCP 

1 1SH  PARAMETER  VALUES.)  RCP 

1117  FORMAT ( 30H  TRY  A LARGER  VALUE  OF  ITERD.)  RCP 

11?0  FORMAT ( 46H  *****  ERROR  *****  THE  VALUE,  RCP 

1 19H  OF  0 FOR  THE  ABOVE)  RCP 

1125  FORMAT (46H  PARAMETERS  LIES  OUTSIDE  THE  ALLOWABLE  RANGE,  , RCP 

1 14HO.OOI  TO  0.35.)  RCP 

1140  FORMAT ( 31H  THE  ERROR  CRITERION  FOR  F IS.E12.2/)  RCP 

1150  FORMAT ( 52H  THF  MAXIMUM  NUMBER  OF  ITERATIONS  ALLOWED  FOR  F IS,  RC» 

1 14/)  RCP 

1160  FORMAT ( 31H  THE  ERROR  CRITERION  FOR  0 IS.E12.2/)  RCP 

1170  FORMAT ( 52H  THE  MAXIMUM  NUMBER  OF  ITERATIONS  ALLOWED  FOR  D IS,  RCP 

1 14/// ) RCP 

1210  FORMAT !5X»1HT,RX,IHP,8X,3HADF«6X, 3HMME , 4X,3HJPA,3X,1HF, BX, 2 HO?,  RCP 

1 8X,2HSF,6X,‘>HMMF(CALC)  ,3X,  5HPHASE/)  RCP 

1220  FORMAT! 3! IX ,F9.4) ,1X,F9.5, 1X.I1 ,2( 1X.F9.6) , IX, IPE9.3,  RCP 

1 21 1X.0PE9.5) ) RCP 

1230  FORMAT ( 53H  *****  ERROR  *****  P MUST  BE  .GE.O.)  RCP 

1240  FORMAT ( 50H  *****  ERROR  *****  IF(JPA.EQ.O)  , RCP 

1 19H  P MUST  BE  .GE.100.)  RCP 

1250  FORMAT ( 50H  *****  ERROR  *****  IF(JPA.EQ.O)  , RCP 

1 17H  p MUST  Bfc  .GF.  T)  RCP 

1260  FORMAT ( 55H  *****  ERROR  *****  T MUST  BE  .GF.-10.IRCP 

1270  FORMAT (55H  *****  ERROR  *****  ADF  MUST  BE  .GE.O.TRCP 

1280  FORMAT ( S4H  *****  ERROR  *****  MF  MUST  BE  .GT.O.)  RCP 

1310  FORMAT (49H  *****  WARNING  *****  POSSIBLY,  RCP 

1 23H  NOT  ENOUGH  ITERATIONS  ) RCP 

1315  FORMAT ( 50H  TO  OBTAIN  SPECIFIED  ACCURACY  FOR  F WITH  THE  ABOVE,  RCP 

1 I’H  PARAMETERS.)  RCP 

1317  FORMAT ( 30H  TRY  A LARGER  VALUE  OF  ITERF.)  RCP 

1320  FORMAT (4HH  *****  ERROR  *****  THERE  IS  NO,  RC® 

I 10H  Rfc  SONANCE ) RCP 

1330  FORMAT! 51H  (DISPLACEMENT  ONE-OUARTER  CYCLE  OUT  OF  PHASE  WITH,  RCP 
1 19H  FORCING  FUNCTION)  ) RCP 

1340  FORMAT ( 25H  FOR  THE  ABOVE  PARAMETFRS)  RCP 

1350  FORMAT ( 51 H *****  WARNING  *****  BECAUSE  OF,  RCP 

1 16H  LARGE  APPARATUS)  RCP 

1360  FORMAT! 49H  DAMPING,  THE  CALCULATED  AMPLITUDE  IS  RELATIVELY,  RCP 

1 15H  INSENSITIVE  TO)  RCP 

1365  FORMAT ( 49H  SPECIMEN  DAMPING.  CONSEQUENTLY,  THE  CALCULATED,  RCP 

l 17H  SPECIMEN  DAMPING)  RCP 

1367  FORMAT! 37H  RATIO  ABOVE  MAY  BE  VERY  INACCURATE.)  RCP 

C SET  INPUT  AND  OUTPUT  OEVICE  COOES  RCP 

ID  * 5 RCP 

00  * 6 RCP 


0056 

0057 

0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 
0081 
0081 
0082 

0083 

0084 

0085 

0086 
00B7 
0088 

0089 

0090 

0091 
009? 

0093 

0094 

0095 
00<>6 

0097 

0098 

0099 

0100 
0101 
010? 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 
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READ{ ID, 10u5 ) EPSF, ITERF, EPSD, I TERD 
IF  IE  PSD.  EQ.  0. ) E PSD  = l.E-A 
IFIITERD.EO.O)  ITERD  = AO 
IFIEPSF.EQ.O.)  EPSF  = l.E-2 
IFIITERF.EC.O)  ITERF  = AO 
WRITE (UDt l 1 AO ) EPSF 
' WRITE tOD, 1150)  ITERF 
WRITE  IOD,  1160  ) EPSO 
WR 1TE  t ODf 1170)  ITERD 
WR  ITE 1 0D»  1 2 10  ) 

10  RE AD  f 1 D , 1000)  T ,P , ADF ,MMF , JPA 
IF  IMMF.EO.O. ) GO  TO  190 
HAL  = 0 

1FIP.GE.0.)  GO  TO  20 
MAL  = MAL  ♦ 1 

IF(MAL.EQ.l)  WRITE! ODt 1 220 ) P,  T, ADF.MMF, JPA 
WRITE (0D,12?0) 

20  IFIJPA.EC.l.OR.P.GE.lOO.)  GO  TO  30 
MAL  = MAL  ♦'  1 

IF(MAL.tQ.l)  WRITEIOD, 1220)  P , T ,A DF ,MMF, JPA 
WR I T£ I OD, 12 AO ) 

30  IF(JPA.eO.l.OR.T.LE.P)  GO  TO  AO 
MAL  = MAL  ♦ 1 

1FIMAL.EQ.1)  WRITEIOD, 1220)  P , T, ADF ,MMF, JPA 
WR  ITE  (OD,  12  SO ) 

AO  1FIT.GE.-10.)  GO  TO  50 
MAL  = MAL  ♦ 1 

IFIMAL.LO.l)  WRIJE 100, 1220)  P , T, ADF ,MMF , JPA 
WRITEIOD, 1260) 

50  IFIADF.GE.O.)  GO  TO  60 
MAL”  = MAC”  VI' 

IFIMAL.EQ.l)  WRITEIOD, 1220)  P , T , ADF ,MMF, JP A 
WRITEIOD, 12  TO ) 

60  IF (MMF.GT.O.)  GO  TO  70 
MAL  = MAL  ♦ 1 

IFIMAL.EQ.l)  WRITEIOD, 1220)  P , T , ADF ,MMF , JP A 
WRITEIOD, 1280) 

70  1FIMAL.NE.0)  GO  TO  10 
KALARM  r 0 
DL  * 0.0001 
DR  = .35 

CALL  R I MI  2 ID, AMPOEL,DELAMP,DL,DR,EPSD,ITERD,IFP ) 
IF(KALARM.Nt.l)  GO  TO  80~ 

MAL  = MAL  ♦ 1 

IFIMAL.EQ.l)  WRITE  IOD,  1220  ) P , T,  ADF  ,MMF  , JP  A 
WRITEIOD, 1320) 

WRITEIOD, 1330) 

WR  I T E I CiD,  1 3 AO  ) 

GO  TO  10 

80  1FIIFR.NE.2)  GO  TO  90 
MAL  = MAL  ♦ 1 

IFIMAL.EQ.l)  WRITEIOO, 1220 ) P, T, ADF.MMF, JPA 
WRITEIOD, 1120) 

WRITEIOO, 1125) 

GO  TO  10 


RCP  0111 
RCP  0112 
RCP  0113 
RCP  011A 
RCP  0115 
RCP  0116 
RCP  0117 
RCP  0118 
RCP  0119 
RCP  0120 
RCP  0121 
RCP  0122 
RCP  0123 
RCP  012A 
RCP  0125 
RCP  0126 
RCP  0127 
RCP  0128 
RCP  0129 
RCP  0130 
RCP  0131 
RCP  0132 
RCP  0133 
RCP  0 13 A 
RCP  0135 
RCP  0136 
RCP  0137 
RCP  0138 
RCP  0139 
RCP  01A0 
RCP  01A1 
RCP  01A2 
RCP  01A3 
RCP  01AA 
RCP  01A5 
RCP  01A6 
RCP  0 1A7 
RCP  01A8 
RCP  01A9 
RCP  0150 
RCP  0151 
RCP  0152 
RCP  0153 
RCP  015A 
RCP  0155 
RCP  0136 
RCP  0157 
RCP  0158 
RCP  0159 
RCP  0160 
RCP  0161 
RCP  0162 
RCP  0163 
RCP  016A 
RCP  0165 
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VO  continue" 

RCP 

0166 

DO  180  1=1,  11 

RCP 

0167 

XR  = (I-D/10. 

RCP 

0168 

OLX  = QL*XR 

RCP 

0169 

PLX  = PL*XR 

RCP 

0170 

CC  = CnSHIQLX )*COSIPLX) 

RCP 

0171 

SC  = SINH(QLX)*C0S1PLX) 

RCP 

0172 

SS  = S1NH(QLX)*SINIPLX) 

RCP 

0173 

CS  = COSH ( OLX )* SIN  C PLX ) 

RCP 

0174 

UPIII  = CI1 )* l-PL*CS+QL*SC>  ♦ CI2)*(- 

PL*SS+OL*CC)  ♦ 

RCP 

0175 

1 C ( 3 )*  I PL*SC  *QL*CS 1 ♦ CI4)*( PL*CC*OL*SS) 

RCP 

0176 

180  CONTINUE 

RCP 

0177 

GAM  = (UPfl>*UP(ll>*4.*!UP<2)*UP(4)+UP<6)*UP(8)+UP(10)) 

RCP 

0178 

1 ♦2.*(UP(3)*UP!5)*UP!7)+UPI9) ))/30. 

RCP 

0179 

SF  = ABS!GAM*F**2/MMF) 

RCP 

0180 

0 = 0*100. 

RCP 

0181 

HR  ITEIOD,  12  20  1 T , P,  ADF  , MMF  , JPA , F , D , SF 

,MMFCAL, PHASE 

RCP 

0182 

1FUER.NE.1  ) GO  TO  110 

RCP 

0183 

WRITE10D, 11 10) 

RCP 

0184 

HR ITEIOD, 11 18) 

RCP 

0185 

WRITE!  OD,  1 1 17  ) 

RCP 

0186 

110  IFIKALARM.NE.2)  GO  TO  120 

RCP 

0187 

HR ITEIOD, 1350) 

RCP 

0188 

WRITEIOD, 1360) 

RCP 

0189 

HR ITE 100, 1365 ) 

RCP 

0190 

WRITE! 00, 1367 ) 

RCP 

0191 

120  IF (KAL  ARM.NE. V)  GO  TO  140 

RCP 

0192 

WRITE! OD, 1310) 

RCP 

0193 

HR  I TE  ( OD,  1315) 

RCP 

0194 

HR  ITE  ! OD,  1317) 

RCP 

0195 

140  CONTINUE 

RCP 

0196 

GO  TO  10 

RCP 

0197 

190  STOP 

RCP 

0198 

END 

RCP 

0199 

Q*********** ************ 

RCP 

0200 

FUNCTION  DFLAMPID) 

RCP 

0201 

c *********************** 

RCP 

0202 

COMMON  /VALS/P.T, ADF, MMF, JPA,AMP,OD,F 

/OLIM/DL.DR 

RCP 

0203 

COMMON  /ALARM/KALARM  /EPSIT /EPS  F,  I TEG  F , EPSO,  ITERD 

RCP 

0204 

EXTERNAL  VFCN 

RCP 

0205 

INtEGER  66 

RCP 

0206 

REAL  MMF 

RCP 

0207 

00  = 6 

RCP 

0208 

DO  = 0 

RCP 

0209 

1FIT.LT.0.)  GO  TO  30 

RCP 

0210 

AA  * — 48*T* P - 7* 1 T*P ) - 1 

RCP 

0211 

Bb  = 48*T*P  ♦ ?0*IT*P)  ♦ 5 

RCP 

0212 

CC  = -4*!T*P)  - 4 

RCP 

0213 

BD2A  = B8/I 2*AA ) 

RCP 

0214 

F = SORT!  12*t  -8 02 A— SORT  I BD2 A**2-CC/AA ) ) ) 

RCP 

0215 

FL  = .3  *F 

RCP 

0216 

FR  = 1.2*F 

RCP 

0217 

GO  TO  60 

RCP 

0218 

30  IF (T.LT.-.l ) GO  TO  40 

RCP 

0219 

FL  = 1.4 

RCP 

0220 

mis  PAGE  IS  BEST  QUALITY  PRACTICABLE 
raOM  COPY  FURNISHED  PO  LDC 
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FR  = 1.4*FL 
GO  TO  60 
40  FL  = 1.75 
FR  = 1.5*FL 
6 0 CONTINUE 

CALL  RTMI  (F  ,V  , VFCN.F L ,FR , EP  SF  , I Tt RF  , I F R » 
IF  (IER.NE.7I  GO  TO  70 
FI  = .99*FR 
FR  = 1.5*FR 

IF  (FL.LE.R.  ) GO  TO  60 
KALAKM  * 1 
GO  TO  120 

70  IF(IFR.EO.l)  KALARM  « 4 
90  CONTINUE 
V = VFCN(F) 

A * AMP 
AM  = MMF/F**2 
OE  LAMP  r A - AM 
IFID.EO.DU  DELL  = DELAMP 
IF 1D.E0.DL)  AML  = A 
IF(O.NF.OR)  GO  TO  120 
PELR  = OELAMP 
AMR  = A 

DIF  = AML  - AMR 

RELOIF  = DIF/AML  ' 

IF (RELDIF.LE.0.20)  KALARM  = 2 
120  CONTINUE 
RE  TURN 
ENO 

C******* *•**•*♦ ******** 

FUNCflON'VFCNlFV  ~ 

C ******* •**•♦♦*****♦** 

DIMENSION  A(4,9),C(4) 

REAL  MMFCAL 

COMMON  /VAL  S/P,T,~S6F  ,MMF,  JP  A,  AMP,D  , FF 
COMMON  /P  AR/PL , OL  ,C 

common  /crit/MMFCal, Phase 

BETA  = SORT (I.*!?. *01**2) 

PL  = F* SORT  I n>ETT*T.  )/2  . I /B  ETA 
Ot  = F*S0RT((BETA-1. 1/2.1/BETA 
SNHO  = SI NH (OL 1 
CSHO  * COSH  (OL  )_ 

SNP  = SIN  (PL  I" 

CSP  * COS (PL) 

CS  * CSHO*SNP  _ 

SC  * SNHO*CSP 
SS  * SNHO*SNP 
CC  « CSMO*CSP 
CSSC  = -PL*C$  *~OL*SC 
SSCC  * -PL*SS  ♦ QL*CC 
sees  = PL*SC  ♦ OL*CS 
CCSS  = PL *CC  ♦ OL*SS 
PF?  = P*F**2 
A ( 1,1)  * T*F**2 
All,2>  * OL  - 2*0*PL 


RCP  0221 
RCP  0222 
RCP  0223 
RCP  0224 
RCP  0225 
RCP  0226 
RCP  0227 
RCP  0228 
RCP  0229 
RCP  0230 
RCP  0231 
RCP  0232 
RCP  0233 
RCP  0234 
RCP  0235 
RCP  0236 
RCP  0237 
RCP  023B 
RCP  0239 
RCP  0240 
RCP  0241 
RCP  0242 
RCP  0243 
RCP  0244 
RCP  0245 
RCP  0246 
RCP  0247 
RCP  0248 
RCP  0249 
RCP  0240 
RCP  0241 
RCP  0252 
RCP  0253 
RCP  0254 
RCP  0255 
RCP  0246 
RCP  0247 
RCP  0258 
RCP  0259 
RCP  0260 
RCP  0261 
RCP  0262 
RCP  0263 
RCP  0264 
RCP  0265 
RCP  0266 
RCP  0267 
RCP  0268 
RCP  0269 
RCP  0270 
RCP  0271 
RCP  0272 
RCP  0273 
RCP  0274 
RCP  0275 
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A 1 1*3)  = - ADF*F**? 

RCP 

0276 

AU, 4)  = i>L  ♦ ?*D*QL 

RCP 

0277 

A ( 2 , 1 ) = -A  ( 1 , 3 ) 

RCP 

0?78 

A 1 2, 2 ) = -A(1,A) 

RCP 

0?79 

AI2,3»  = A( 1, 1 | 

RCP 

0280 

A ( 2 , A 1 = A ( 1,2) 

RCP 

0281 

A(3,l)  = CSSC  - 2.*D*SCCS  - PF2*CC 

RCP 

028? 

AO, 2)  = SSCC  - 2.*D*CCSS  - PF2*SC 

RCP 

0283 

AO,  3)  * SCCS  ♦ ?.*D*CSSC  - PF2*SS 

RCP 

028  A 

AI3.A)  = CCSS  ♦ 2 .*0* SSCC  - PF2*CS 

RCP 

0288 

A C A , 1 ) = -A  ( 3 , 3 1 

RCP 

0286 

AI-,2)  = — A 1 3 , A ) 

RCP 

0287 

A(A,3)  = A(3, 1) 

RCP 

0288 

AIA.A)  = AO,?) 

RCP 

0289 

Cll)  = 0 

RCP 

0290 

CC?>  = -1. 

RCP 

0291 

C C 3)  = 0 

RCP 

0?9? 

CO)  * 0 

RCP 

0293 

EPS  = l.E-5 

RCP 

029A 

NN  = A 

RCP 

0298 

CALL  GELG(C, A , NN, 1 , E PS , I ER ,NN* NN , NN ) 

RCP 

0296 

UO  = C ( 1 ) 

RCP 

0 297 

VO  * C (3) 

RCP 

0?°8 

AMPLO  = SORT  I UO**?  ♦ VO**?) 

RCP 

0299 

PHASEO  = A T AN2  C VO  , UO  ) 

RCP 

0300 

IF (PHASEO.LT.O. ) PHASEO  = PHASEO  ♦ 6.283185 

RCP 

0301 

UL  * C(1)*CC  ♦ C ( 2 ) *SC  ♦ C ( 3) *S S ♦ C(A)*CS 

RCP 

030? 

VL  = C ( 3 ) *CC  ♦ C t A ) *SC  - C(1)*SS  - Ct?)*CS 

RCP 

0303 

AMPLL  * SQR  T (UL **?  ♦ VL**?) 

RCP 

030A 

PHASEL  = ATAN2IVL.UL  ) 

RCP 

0308 

IF(JPA.EO.O)  VFCN  * PHASEO  - 3.1A1593 

RCP 

0306 

IF(JPA.NE.O)  VFCN  * PHASEL 

RCP 

0307 

IF(JPA.FQ.O)  AMP  = AMPLO 

RCP 

0308 

IF(JPA.NE.O)  AMP  * AMPLL 

RCP 

0309 

MMFCAL  = AMP*F**2 

RCP 

0310 

PHASE  = VFCN 

RCP 

0311 

RETURN 

RCP 

0312 

END 

RCP 

0313 

RCP 

031 A 

SUBROUT INE  RTMI2(X,F,FCT,XLI.XRI,EPS,IFND,IER) 

RCP 

0315 

COMMON  /ALARM/KALARH 

RCP 

RCP 

0?16 

0317 

I E R=0 

RCP 

0318 

XL=XLI 

RCP 

0319 

XR  = XR I 

RCP 

0320 

X*  XL 

RCP 

0321 

TOL*  X 

RCP 

0322 

F*FCT( TOL ) 

RCP 

0323 

IF (KALARM.FO. 1 ) RETURN 

RCP 

032A 

IF ( F ) l , 16, 1 

RCP 

0325 

1 fl=f 

RCP 

0326 

X*  XR 

RCP 

0327 

TOL*  X 

RCP 

0328 

F*FCT«  TOL  ) 

RCP 

0329 

IF (KALARM.EO, 1 ) RETURN 

RCP 

0330 
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I F ( F ) 2 f 16,2 

RCP 

0331 

2 

FR=F 

RCP 

0332 

If <S1GN(1.,FL1*SI6N|1.,FR> )25,3 .25 

RCP 

0333 

3 

1*0 

RCP 

0334 

TOLF=100.*EPS 

RCP 

0335 

* 

1 = 1*1 

RCP 

0336 

DO  13  K = 1 » I END 

RCP 

0337 

X= .5* 1 XL+XR ) 

RCP 

0338 

TOL*X 

RCP 

0339 

F=FCTI TOL) 

RCP 

0 340 

IF  (KALARM.EQ.  1 ) RETURN 

RCP 

0341 

IF(F)5, 16,5 

RCP 

0342 

5TF  (SIG  NTIT,T7T5TG>i7  IT,  FR  )) 7 , 6 , 7 

RCP 

0343 

6 

TOL=XL 

RCP 

0344 

XL=XR 

«CP 

0345 

XR=TDL 

RCP 

034b 

TOL=FL 

RCP 

0347 

FL=FR 

RCP 

0348 

FR=TOL 

RCP 

0349 

7 

TOL=F-FL 

RCP 

0350 

A=F*TOL 

RCP 

0351 

A=A*  A 

RCP 

0352 

Tf  (A-FR*( FR-FL) >8,9,9 

RCP 

0353 

6 

IFII-IFND)17,17,9 

RCP 

0354 

9 

XR=X 

RCP 

0355 

FR  =F 

RCP 

0356 

TOL*f PS 

RCP 

0357 

A=ABS(XR) 

RCP 

0358 

I F ( A-l . ) 1 1 , 1 1 , 1 0 

RCP 

0349 

10 

TOL=TOL*A 

RCP 

0360 

' ii 

I F ( A 6 S ( XR— XLT— T OL  1 1 2 , 12,13 

RCP 

0361 

12 

IF l ABS  ( FR— FL)— TOLF) 14, 14, 13 

RCP 

0362 

13 

CONTINUE 

RCP 

0363 

1ER=  1 

RCP 

0364 

14 

IF  (ABS«FR>-ABS(FU)  16,16,16 

RCP 

0366 

15 

X = XL 

RCP 

0366 

F = FL 

RCP 

0367 

lb 

RE  TURN 

RCP 

0368 

17 

A=FR-F 

RCP 

0369 

DX  = (X-XLI*FL*  C1.*F*C A-TOL)/ (A*( FR-FL) I) /TOL 

RCP 

0370 

XN  = X 

RCP 

0371 

FM=F 

RCP 

0372 

X*XL— DX 

RCP 

0373 

TOL*  X 

RCP 

0374 

F=FCT( TOL ) 

RCP 

0375 

IFIKALARM.EO.il  RETURN 

RCP 

0376 

IF (F 116,16, 18 

RCP 

0377 

18 

TOL=EPS 

RCP 

0378 

A* AB  S ( X I 

RCP 

0379 

IFIA-l. 120,20,10 

RCP 

0380 

19 

TOL=TOL*A 

RCP 

0381 

20 

IF ( ABS  ( DX I— TOL ) 21 ,21,22 

RCP 

038? 

21 

IF(ABS(F)-T0LF)16,16,22 

RCP 

0383 

22 

IF ( SIGN ( 1., F ) *SIGN(1.,FL II 26,23,24 

RCP 

0384 

23 

XR=X 

RCP 

0385 
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FR=F 
GO  TO  4 

24  XL  =X 
FL  = F 
XR=XM 

_FR  =FM 
GO  TO  4 

25  IER=2 
RETURN 
END 


RCP  0386 
RCP  0387 
RCP  0388 
RCP  0389 
RCP  0390 
RCP  0391 
RCP  0392 
RCP  0393 
RCP  0304 
RCP  0395 
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r 


F ! LF  ; G^LGOTv)!  SUBRT  A 


CONVFRS AT  1 ONAL  MONITOR  SYSTEM 


subroutine  gels  ir,a,m,n,eps,  ier.mm.mnj 

RCP  0396 
RCP  0397 

r *******  ************  ******************  ******* 

DIMENSION  A(MM),  R(MN) 

nrP  039R 
RCP  0399 

IF<M|23,23,1 
L l5°=o 

RCP  0400 
RCP  0401 

PI V=0. 

RCP  0402 

RCP  0403 

RCP  0404 

DO  3 L=1,mm 

RCP  0405 

TB=A8$1ACL)I 

RCP  0406 

IF(TR-PIVI3.3.2 

RCP  0407 

2 PI V=TR 

RCP  040R 

I = L 

RCP  0409 

3 CONTINUE 

RCP  0410 

TOL  = Fp  S*p  I V 

RCP  0411 

LS  T = 1 

RCP  0412 

DU  IT  K=1.M 

RCP  041,3 

IF(PIVI23,23,4 

RCP  0414 

4 I F ( IER)7,5,7 

RCP  0415 

5 IF(PIV-TOL 16,6,7 

RCP  0416 

6 ICR=K-1 

DTP  0417 

7 P1VT=1./A(  n 

RCP  0418 

j=<  r-t i/m 

RCP  0419 

!=  T-J*M-K 

RCP  0420 

J= J+J-K 

RCP  0421 

OP  9 L=K,NM,M 

RCP  0422 

LL  =L  ♦! 

RCP  0423 

Tn=PIVI*R( III 

PC P 0424 

o(LLI=PCU 

°CP  0425 

fl  R/L»=T8 

RCP  0426 

I F ( K-M )9,18,1R 

PfP  04?7 

9 LEND=L  ST  «-m-k 

RCP  0428 

Ic( JI12, 12,10 

°CP  0429 

10  II=J*M 

PC®  0430 

DO  11  L=LST,LFNO 

RCP  0431 

T R = A ( L ) 

RCP  0432 

LL  =L  ♦ I I 

RCP  0433 

At  11  = A (LL 1 

RCP  0434 

11  A(LLI*TR 

RCP  0435 

12  OP  13  L = L S T , MM , M 

RCP  0436 

LL  = L ♦ ! 

*>CP  043^ 

TR  = PtV  I*.M  LL  t 

RCP  0438 

A(LLI=A(L» 

RCP  0439 

13  a(li=tb 

RCP  0440 

A(t?n  =j 

#rp  0441 

p I v =Y5 . 

RCP  0442 

LST=LST*1 

RCP  0443 

J*0 

SCP  0444 

OP  16  II=LST,LFND 

RCP  0445 

pTvt=-a{ n j 

■Rfp  04  46 

I$T=II*m 

RCP  044'’ 

j*  j*T 

RCP  044S 

00  15  L=1ST,MM,M 

RCP  0449 

LL=L-J  PC P 0450 


1 
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THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
COPY  FUptilSHEB  TO  BD,C  - 


FILE:  GELGRTMI  SUP»t  A CON VERS AT IQNAL 

monitor  system 

MLi  = MLi*pTvT*ML11 

PCP 

0461 

TB=ABS 1 A(L  1 1 

RCP 

0452 

If IT8-PI V) 15,15, 14 

RC«> 

0453 

14  P I V = T R 

RCP 

0454 

1=1 

RCP 

0455 

IS  CONTINUE 

RCP 

0456 

00  16  L=K,NM,m 

PC° 

0457 

LL=L*J 

PC  p 

0453 

16  R(lLI=8(LL)+PIVI*5(Ll 

RCP 

0454 

17  LST*LST*M 

RCP 

0460 

i a tttm-i 12?, 22,16 

rTP 

7461 

19  IST=MM*M 

RCP 

0462 

LSt=M*t 

6£P~ 

046? 

00  21  1=2, M 

RCP 

0464 

1 1 =LST- l 

RCP 

0465 

I ST= IST-LST 

RCP 

0466 

L=  I ST-M  “ 

RCP 

046'' 

L=AtLI».5 

RCP 

0468 

DO  21  J= I I , NM , M 

RCP 

0469 

TR=R(J 1 

RCP 

0470 

LL  = J 

RCP 

0471 

0^  20  K= I S T , MM , M 

RCP 

0472 

LL  =L  L* l 

RCP 

0473 

20  TB=T8-A(K)*R<LLI 

RCP 

0474 

K=  J + L 

RCP 

0475 

RC J)=R (K  » 

RCP 

0476 

FT  RTkT=tr 

RCP 

04T? 

22  RFTllRN 

RCP 

0478 

7?  Te*  =-l 

RC  ° 

0479 

RETURN 

RCP 

0430 

C*IT5 

HTTP" 

f ************************************************** 

RCP 

0482 

5 URROUTINF  RT  ,vj  ( X,  F,  F C T , xrrTXPTTFPS , IcNO,IER| 

RC° 

0433 

C *************************************  ****** ******* 

RCP 

0484 

IE  6 = 0 

“6FP 

04?  5 

XL  = XL  1 

RCP 

0436 

X°  = XR 1 

RCP 

74  3? 

X=  XL 

RCP 

0483 

T5l*x 

IFF6 

7479 

F>FCT( TOLI 

RCP 

0490 

THTTl ,16,1 

P76" 

7491 

1 FL  = F 

RCR 

0492 

x = y o 

RCP 

0493 

TOL  = X 

RCP 

0494 

F=Fctit0l7 

B7P 

0495 

IF(f)2,16,2 

RCP 

0496 

T fp*f 

RCP 

7477" 

I F ( S IGN ( la ,FLI*SIGN1 1..FR1 125,3,25 

RCP 

0498 

3 1 = 0 

RCP 

w>L%:L-M 

T0LF»100.*EPS 

RCP 

0600 

4 T*  T*T 

HTTP" 

'0501 

00  13  K*  1 , IEND 

RCP 

0502 

"11=4  5 » 1 XL  ♦ X R 1 

TOL  =X 

RCP 

0504 

F = FC"H  TULJ 

RCP 

05U5 

THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
roOM  COPY  BURNISHED  TO  liDC 


FJL F:  GCLGRTHI  SU8RT 


C ON VERS AT  I ONAL  MONITOR  SYSTFH 


IF<F»5,16,5 

PCP 

0506 

5 

PCP 

05  OT 

6 

TOL=XL 

Rf.P 

0508 

XL  = XR 

RCP 

0509 

XR  = TOL 

RCP 

0510 

TOL=FL 

PCP 

0511 

fl*fr 

PCP 

0512 

FR  = TQL 

PCP 

0513 

7 

tol=f-fl 

PCP 

0514 

A=F*T3L 

RCP 

0515 

A = A*A 

PCP 

0516 

IFIA-F9*(FR-FLI  >8.9,9 

PCP 

05  1 7 

8 

I F ( I-IENDI17.17.9 

RCP 

0518 

9 

XP  =x 

RCP 

0519 

FR  = F 

RCP 

0520 

T0L=cPS 

RCP 

0521 

A=  ABS( XRI 

PCP 

0522 

I F(  A-l  • ) 11.11,  10 

RCP 

0523 

10 

TOL*TOL*A 

RCP 

0526 

11 

I F ( ABS  C XP— XL l-TOL) 12,12,13 

RCP 

0525 

12 

IF ( ABS1FR-FL l-TOLF 116,14,13 

RCP 

0K26 

13 

CONTINUE 

RCP 

0527 

IF  R=  1 

RCP 

0528 

1* 

IF (ABSCFRI -ABS fFLI 116,16,15 

PCP 

0529 

15 

X=  XL 

PCP 

0530 

f*fl 

PCP 

0531 

16 

R 5 TURN 

RCP 

0532 

17 

A=  F P-F 

RCP 

0533 

DX=( X-XL  >*FL*( l.*F*| A-TOLI/I A*IFR-FD) l/TOl 

PCP 

0534 

XM  = X 

RCP 

0535 

RCP 

0536 

X=  XL-OX 

RCP 

0537 

TOL  *X 

PCP 

0538 

F=  FC T(  TOL  1 

RCP 

0539 

IF  C FI18, 16,18 

PCP 

0540 

18 

TOL=FPS 

RCP 

0541 

A=A8S(  XI 

RCP 

0542 

1FCA-1. 120,20, 19 

RCP 

0543 

19 

tol=tol*a 

RCP 

0544 

20 

IFIABSIDXI-T OLI21,2l,22 

RCP 

0545 

21 

I c ( ABSI FI-TOLF 116, 16,22 

6£p 

Mil 

22 

Ic 1 5 IGNI 1 . , F 1 +S I GN ( l.,FL  1 1 24,23,  24 

RCP 

054T 

23 

XR  = X 

RCP 

0548 

FR  *F 

RCP 

0569 

GO  TO  4 

RCP 

0550 

26 

XL  = X 

RCP 

0551 

FL  = F 

SfP 

[■VIM 

XP  = XM 

RCP 

0553 

FR  * FM 

5CP“ 

GO  TO  4 

RCP 

0555 

rr 

IFR*2 

0556 

return 

RCP 

0557 

INT5 RCP  055S 
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RESONANT  COLUMN  TEST 
(Data  Shaet) 


Technician  Computed  By  Checked  By 


RESONANT  COLUMN  TEST 
(Intermediate  Calculations  Sheet) 


Technician  Computed  By  Checked  By 


RESONANT  COLUMN  TEST 
(Final  Calculations  Sheet) 


Technician  Computed  By 


RESONANT  COLUMN  TEST 
(Amplitude  Decay  Data  Sheet) 


PLATE  4 


Technician  Computed  By  Checked  By 


In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Drnevich,  Vincent  Paul 

Resonant  column  test  / by  Vincent  P.  Drnevich,  Soil  Dynamics 
Instruments,  Inc.,  Lexington,  Kentucky.  Vicksburg,  Miss.  : 

U.  S.  Waterways  Experiment  Station  ; Springfield,  Va.  : avail- 
able from  National  Technical  Information  Service,  1978. 

63,  c7^  p.,  4 leaves  of  plates  : ill.  ; 27  cm.  (Miscellaneous 
paper  - U . S.  Armv  Engineer  Waterwavs  Experiment  Station  ; 
S-78-6) 

Prepared  for  Office,  Chief  of  Engineers,  U.  S.  Army,  Wash- 
ington, D.  C.,  under  Contract  No.  DACW39-77-M- 1687 . 

Includes  bibliographical  references. 

1.  Computer  programs.  2.  Damping.  3.  Resonant  column  tests. 

4.  Shear  properties.  5.  Soil  test  specimens.  6.  Vibrations. 

1.  Soil  Dynamics*  Instruments.  Inc.  II.  United  States.  Army. 
Corps  of  Engineers.  111.  Series:  United  States.  Waterways 
Experiment  Station,  Vicksburg,  Miss.  Miscellaneous  paper  ; 
S-78-6. 

TA7 . W34m  no. S-78-6 


